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THE  GENERALIZED  SIMPLE  WAVE1 
by 

D.  Naylor2 

Introduction.  An  important  exact  solution  in  the  aero- 
dynamics 0f  steady,  perfect  and  irrotational  compressible  flow 
is  the  Prandtl-Meyer  expansion  wave.  This  solution  was  originally 
obtained  by  seeking  directly  a solution  of  the  equations  of  com- 
pressible flow  referred  to  polar  coordinates  in  which  the  radial 
distance  is  of  no  consequence  to  the  density,  pressure  and 
magnitude  and  direction  of  the  velocity.  It  therefore  had  the 
property  that  the  density,  pressure  and  velocity  vector  were 
constant  along  radii  through  a fixed  point,  but  this  restriction 
can  be  removed  and  the  solution  is  capable  of  generalization  to 
the  case  for  which  conditions  are  constant  along  a series  of 
straight  lines  not  necessarily  concurrent.  In  the  extended  case 
the  variables  u,  v,  p,  p are  all  functions  of  just  one  independent 
variable  and  relations  of  the  type  3 (u ,v)/9 (x,y ) = 0 hold  every- 
where. The  vectors  Vu,  Vv,  Vp,  vp  are  all  parallel  and  each  of 
the  variables  may  be  considered  as  reducing  to  a function  of  a 
function  X(x,y)  which  may  remain  unspecified  or  chosen  to  reduce 
to  q,  or  8 (say),  or  any  combination  of  these  quantities,  as 
convenience  requires.  Such  flows  for  which  the  velocity  vector 

371  The  results  presented  in  this  paper  were  obtained  in  the 

course  of  research  sponsored  by  the  Office  of  Naval  Research 
under  Contract  NVonr-35801  with  Brown  University. 

» Research  Associate,  Division  of  Applied  Mathematics,  Brown 
University. 
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2.  reduces  to  a (vector)  function  of  just  one  independent  variable 
X are  now  commonly  referred  to  as  "simple  waves".  It  is  known 
that  steady  irrotational  flows  defined  in  this  way  have  the 
property  that  the  pressure  and  density,  as  well  as  the  velocity 
q,  reduce  to  functions  of  X, 

More  generally  it  is  desirable  to  extend  the  concept  of 
the  simple  wave,  for  which  q reduces  to  a function  of  X only, 
but  to  say  nothing  about  the  remaining  variables  p and  p , to 
flows  which  are  not  necessarily  steady  and  irrotational  or  for 
which  conditions  of  velocity  are  not  necessarily  constant  along 

straight  lines  or  that  the  isovels  intersect  in  a fixed  point, 

* 

The  first  step  in  this  direction  was  taken  by  Prin  [1]  who  con- 
sidered steady  rotational  plane  flow  by  adapting  the  Prandtl- 
Meyer  method  thereby  presupposing  concurrent  isovels.  Martin, 
in  a series  of  lengthy  investigations,  [2,  3,  h)  succeeded  in 
dealing  with  the  more  general  case,  without  making  any  assumptions 
about  the  isovels,  by  extracting  the  simple  wave  from  a general 
theory  in  which  the  pressure  p and  the  stream  density  function 
were  selected  as  fundamental  variables.  It  appeared  that  the 
isovels  constituted  a system  of  straight  lines,  parallel  or  con- 
current, Ghaffari  [ 5 1 applied  the  original  steady  irrotational 
simple  wave  theory  in  terms  of  q as  the  fundamental  variable  with 
the  aid  of  the  Legendre  contact  transformation  and  Giese  [6] 
later  used  this  same  procedure  in  terms  of  X to  discuss  essen- 
tially geometrical  properties  of  the  hodograph  for  one  precisely 

* Numbers  in  square  brackets" refer  to  the  Bibliography  at  the 
end  of  the  paper. 
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similar  phenomena  including,  in  addition  to  the  simple  wave 
q(X),  the  "double"  wave  q(X,p). 


The  One  Parameter  Fields.  Only  when  the  flow  is  unsteady 

is  the  velocity  a function  of  two  independent  variables  in  which 

case  the  problem  is  mathematically  a double  wave,  although  the 

flow  physically  is  instantaneously  a simple  wave  varying  with  the 

time.  The  one  parameter  velocity  field  is  a special  application 

of  a more  general  method  in  which  any  property  of  the  flow  may 

be  regarded  as  being  reducible  to  a function  of  only  one  parameter. 

In  unsteady  flow  the  method  may  be  employed  by  focusing  attention 

on  velocity  potentials  of  the  type  cp(X)  or  more  generally  q>(X}t). 

In  linearized  flow  the  method  yields  supersonic  potentials  when 
2 p p P p 

X'  = x - £ (y  + z ).  The  reduced  potential  of  oscillatory  flow 
is  of  this  variety, 


f (x,y,z) 


cos  n0X 

£ 


2 P 2 

with  the  standard  notation  B = M - 1,  nD  = u/aB  » The  curves 
X(xyt)  = constant,  are  the  characteristic  curves  which  in  the 
generalized  simple  wave  reduce  in  all  cases  to  a system  of  straight 
lines  (not  necessarily  coincident  with  one  family  of  Mach  lines 
of  supersonic  flow)  except  in  the  rotational  unsteady  compressible 
case. 


The  purpose  of  this  report  is  to  show  how  a variation  of 
the  hodograph  method,  based  on  the  application  of  the  Legendre 
contact  transformation  in  terms  of  the  parameter  X which  remains 
unspecified,  may  be  conveniently  employed  to  rapidly  construct 
rotational  anc.  unsteady  generalizations  of  the  simple  wave  q(X) 
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when  X = X(x,y,t)  in  general.  Finally  the  more  general  unsteady 
wave  for  which  ^ = q(X;t)  5 X = X(xyt)  is  considered  in  Section  IV. 

General  Equations  of  the  Flow.  The  general  equations 
of  motion,  continuity  and  adiabatic  transformation, 

("at  + i * + p"-1  vp  = 0 


f!  + v • (pa)  - ° 

(■§£  + q • v)  (pp~Y)  " 0 


may  be  transformed  to  give 


2 _2_ 

a , Y“l»  , Y-l 

(aT  , + v * a ) = 0 


at 


(1.1) 


where 


3d 

+ VH  - TVS  = q W 

at  -A  “ 


S = Cy  log  (pp  Y) 


(1.2) 


W = va2 

1/2 

and  a = (yp/p)  defines  the  acoustic  speed  and  T = p/(pR)  the 

p 

absolute  temperature.  Since  T - a /yR  then  Crocco's  equation 

(1,2)  may  be  used  to  yield  another  equation  in  addition  to  (1.1) 

2 

involving  just  a and  q 

M3'2  (is  + a • v)a)  = 0 a.3) 

or 

V.-la"2  (q,  H-  qVl  - q.  W)|  =0.  (1.U-) 

At  —T, 
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Steady  Rotational  Flows;  A Simplified  Theory.  Attention 
being  confined  to  iso-energetic  flows  for  which  H = constant, 


the  general  equations  for  the  velocity  reduce  to 

2 

V (q  aY  1)  = 0 


where 


VA(a"2  q AW)  = 0 

2 ? 2 
q . a"1  . _ c^ 

T + 7^1  “ “2' 


(1.5) 

(1.6) 

(1.7) 


serves  to  provide  the  acoustic  speed  directly  in  terms  of  the 
velocity  and  the  iso-energetic  constant  c (actually  the  greatest 
possible  speed).  The  fundamental  assumption  of  the  existence 
of  a flow  of  the  simple  wave  pattern  is  most  easily  exploited 
by  using,  not  the  . rdinary  stream  function,  but  Crocco’s  stream 
function  \|>  given  by  2 


V = Ka  q a 


Y-l 


(1.8  i 


the  necessary  condition  for  the  existence  of  which  is  implied  by 
(1,5)}  and  then  define  the  associated  function  ®(x,y)  by  the 
Legendre  contact  transformation 


I 


i|)  = ?'  + r • 2. 

_,2 


(1.9) 


with  the  notation  £ = q a 1 ~ = Vty  which  in  the  following 
analysis  is  convenient  because  £ rather  than  q,  is  the  fundamental 
field  vector. 

The  transformations  employed  are  two-fold i 
1,  The  current  independent  variables  (x,y)  are  replaced  by 
<M)t 
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2,  The  dependent  function  ij>(x,y)  is  itself  replaced  by 

its  Legendre  transform,  but  is  also  retained  as  a fund- 
amental independent  variable  where  necessary. 

This  procedure  reduces  the  general  problem  of  the 
rotational  simple  wave,  originally  that  of  two  non-linear  partial 
differential  equations  for  q. ( X) ; w(x,y)  to  the  study  of  an 
ordinary  second  order  differential  equation  for  q(X)  and  a 
linear  ordinary  second  order  differential  equation  for  T’CX). 

This  reduction  is  rapidly  achieved  because  £ and  ? both  reduce 
to  functions  of  X only,  they  are  not  general  functions  of  the 
new  variables  (X,\j>),  as  for  example,  is  w(X,i|j),  The  transforma- 
tion (Xj'.y) — > (X,40  is  permissible  because 


a(XTilO 

9(x,y) 


* 0 


which  independence  is  immediate  because 


VA  - - axA  V*.  - o 

in  virtue  of  = 0. 

It  follows  that  VX.  - 0. 

A 

that  is,  £ VX  = 0 

only  when  £ £ ~ 0 

A X 

which  is  the  excluded  degenerate  case,  9^  ~ 0.  of  a parallel 
(shear)  flow.  The  following  outlines  and  illustrates  the  general 
method,  the  principles  of  which  are  subsequently  adapted  to 
study  unsteady  flows  in  the  sections  which  follow*  Full  details 
and  results  of  the  rotational  one-parameter  fields  will  be  found 
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in  the  author’s  memoir  [7]  where  a slightly  different  treatment 
is  employed.  It  is  emphasized  that  the  adopted  procedure  is  not 
restricted  to  flows  of  the  type  here  discussed,  but  with  modifi- 
cation is  capable  of  application  to  a great  diversity  of  fluid 
flow  phenomena  of  the  simple  wave  pattern,  (e,g.  unsteady, 
rotational  and  viscous  incompressible  motion  - see  another 
report),  the  criterion  being  the  existence  of  a suitable  field 
gradient  vector  depending  on  X only.  It  is  at  once  clear  from 
(1.7),  and  (1.8)  that  a , and  2,  reduce  to  functions  of  X only, 
being  independent  of  i|j,  therefore  it  is  possible  to  write  the 
total  differential  of  (1.9)  in  the  form 

d<P  + (r  * &x)dX  = 0 

so  that  ? = 'P(X)  and 

+ r ’ 2X  = 0 (1.10) 

from  which  by  further  differentiation  is  obtained 


(fxx  + - ■ ^xx)vx  + ^x  = 


= 0 


(1.11) 


The  following  representations  are  employed, 

2Y 


| q | = D(X)  a 
“>A  * 


r-i 


(1.12) 


where 


+ F(X)0  = K(X)& 

A.A  \ 

ah)(X)  = (a^  - q^)q^  + a^q^©^ 

F(X)  = ■ aT  (los  V 

p 3q2  2 


(1.13) 

(1.14) 

(1.15) 


(1.16) 
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Also  an  expression  for  the  vorticity  which  follows  from  (1.6), 
and  (1.8) 

- ^r 

Vtb  A V(w  a Y-l)  = 0 

A 2r 

or  TTT 

w = G(t|j)a^  (1.17) 

in  terms  of  some  function  G(tJj). 

Since  w q = V\  q.  then  (1,11),  (1.12), and  (1.17) 

gives 

0i>xx  + r G (^)  = D(\) 

finally  substitution  of  the  scalar  produce  £ • £ . from  (1.13) 
with  the  aicl  of  (1.9),  and  (1.10)  gives 

[3>x>  + Fip^  + K( ^ - ?)  ] G (xl>)  = D(  X)  (1.18) 

Since  (X,-|)  are  independent,  a field  equation  of  the  type  (1,18) 
can  only  be  true  (since  F,  K,  ?,  D are  independent  of  i|>)  if 
either 

1.  G (ij> ) - A | K(X)  = 0;  9>xx  + F(X)?X  = 2 1^2. 

2,  G(\jj)  = A/\Jj | AK(X)  = D(\)  5 ?>xx  + F(\)3>x  - K(\)'P  = 0 

where  A is  a constant  determining  the  vorticity  distribution  in 

alternative  forms  2y  2y 

Y-l  v-T 

A a • (A/^)  a' 

in  terms  of  Croccots  stream  function  if  formulae  expressed  in 
terms  of  the  normal  stream  function  are  required  then  appropriate 
conversion  relations  must  be  employed. 
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Properties  of  the  X-lines.  The  first  possibility 
K(\)  = o implies,  from  (1.13), 

fiu  + P(X)%  “ 0 

which  may  be  integrated  to  show  that  £x  is  parallel  to  a fixed 
vector*  Therefore  the  X-lines  (1,10)  constitute  a system  of 
parallel  lines.  These  are  rhe  flows  which  appear  to  elude  treat- 
ment by  other  general  methods  and  for  which  a direct  approach 
is  simpler* 

The  second  possibility,  for  which 

tpx\  + - K(X)T  = 0 (1,19) 

+ F(X)\  " K(X)&  = ° (1.20) 


implies , 


eliminating  the  functions 


I 


<P 


XX 


u 

ux 

uxx 


F(X),  K(X) , 
V ' 

I 

vx  l = ° 
vxx  ! 


where  (U,V)  are  the  components  of  C3(X),  This  means  however  that 
the  three  non-parallel  lines 


<P-rr*£  = 0;  «PX  + v * £x  = 0, 


?XX+* 


0 


are  concurrent  and  implies  the  existence  of  a vector  r 0(M  such 
that 


171  4- 


^o 


i n - 
SL  ~ 


A. 


- v»  -A 

*t>'  *X 


— A • 

— W r> 


XX 


+ r 
“O 


n 

*XX 


identically  in  X, 
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Differentiating  the  first  one  and  using  the  second  gives 

<Eo>x  • 2 = 0 

whilst  differentiating  the  second  and  using  the  third  yields 


(£0>x  • = o 

so  that  since  the  vectors  £>,  £ are  neither  zero  nor  parallel 
(except  in  degenerate  cases). 

IT  “ 

Therefore  r0  is  a fixed  vector,  independent  of  X,  and  the  three 
lines  in  question  all  pass  through  the  same  fixed  point  for  all 
values  of  X;  in  particular  the  parameter  lines  (1,10)  are  con- 
current, These  results  shovr  that  in  no  case  is  a non-degenerate 
envelope  of  the  X-lines  possible,  and  also  that  in  the  case  of 
concurrent  X-lines,  that,  by  selecting  as  the  origin  the  point 
through  which  they  all  pass,  'P(X)  = 0,  In  other  words,  the 
Legendre  transform  of  Crocco’s  stream  function  may  be  taken  to 
be  identically  zero  in  the  theory  of  a rotational  steady  simple 
wave,  except  when  conditions  are  constant  along  a series  of 
straight  lines.  This  is  in  direct  contrast  to  the  classical 
simple  wave  theory  in  steady  potential  flow  for  'which  an  arbitrary 
convex  X -envelope  is  possible  corresponding  to  arbitrary  ’-PCX). 
The  preceding  method  is  still  valid,  however,  as  far  as  equation 
(1.18),  which  in  the  simpler  case  of  G(^»)  * Q simply  reduces 
to  D(X)  = 0 so  that  ’P(X)  may  be  arbitrarily  chosen  (subject  to 
the  condition  that  the  X-envelope  is  a simple  convex  curve).  By 
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(1.1*0  the  equation  for  the  velocity  D(\)  = 0 will  yield  by 
integration  the  standard  hodograph  epicycloid  0 = 0(q).  More 
directly  (1,5)  gives  2 

V X • (q  aY"1)x  = 0 

and  y.  q = 0 gives 

V\A  q = 0 
A 

therefore  the  velocity  relation  for  a classical  simple  wave  is 

2 

X 

* (q  a'  0 

which  on  expansion  yields  D(X)  = 0, 

2,  The  Simple  Wave  in  Unsteady  Flow,  When  unsteady  flow 
is  considered  a stream  function  never  exists  in  compressible  flow 
and  it  is  necessary  to  confine  attention  to  irrotational  motion 
in  order  to  adapt  the  preceding  method  by  replacing  the  stream 
function  of  Crocco  y by  the  velocity  potential  function  cp(x,y,t). 
The  equations  of  motion,  continuity,  irrotational  motion  and 
isentropic  flow  are 

pt  + v • (pq)  = o 

q.  + (q  • v)q  + 1 vp  ~ o 

~t  p 

VA  - = ° 

„Y 

pp  = constant 


may  be  used  to  yield  the  general  equations  of  unsteady  potential 
flow  in  the  form 


;«w 


f 


*11-111 


" (ab)  + V * (q  ab)  = 0 

2 2 «_  „2 
SL  + a.  + = £_ 

2 y I 9 t 2 
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(2,1) 

(2,2) 


where  c is  a constant,  b = 2/y-l,  q = Vcp , and  a = (^2.)b//2  is  the 
acoustic  speed, 

Introd\icing  <I>,  the  Legendre  transform  of  9, 


then 


(p  = <J>  + r • Vcp 


d$>  + r • dq  =0 


(2.3) 


The  simplest  generalization  of  the  simple  wave  will  be  one  for 
which  the  velocity  vector  q reduces  to  a function  of  just  one 
independent  parameter  \ incorporating  the  time,  X =X(r  ; t)» 
Therefore  q = q(>0  and  dq  = q^d\  so  that 

d$  + r 4£^dX  = 0 

V$  + (r  • q )7\  = 0 

""" 

therefore  $ = <±>(X,t)  and 

+ - ’ = 0 (2.1!-) 

the  function  $ reducing  to  a function  of  (\,t)  only,  involving 
t independently  as  well  as  through  \(r  ; t). 

The  surfaces  of  constant  \ , that  is  of  constant  q,  jt 
constitute  a system  of  planes  which  vary  with  the  time.  In  plane 
flow  the  velocity  is  constant  along  each  of  a family  of  straight 
lines  varying  with  the  time  and  normal  to  q^,  The  angle  of 
intersection  T(\)  with  the  streamlines  being  given  by 


cot  T = (q8^)/q^ 


(2.5) 


VP5 


. •*  1 


V 
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are  equivalent.  When  operating  on  q(X)  this  operator  reduces 
to  3/3\;  but  the  general  form  must  be  retained  for  functions 

f ..  n 

like  9,  9 and  a^  = yp/p»  since  they  will  in  general  be  functions 
of  9 as  well  as  (X,t)« 

Actually  in  the  one  parameter  fields  here  considered 
only  9 does  not  in  general  reduce  to  a function  of  (X,t)  only* 
The  equations  (2*1)  and  (2,2)  for  the  velocity  now 
become  2 


(a^q  - qq^q)  * VX  = 2qn  X + 9 

2 2 2 
SL  + + <S?  = ° 

2 + Y-l  t ^ 

• *• 

and  the  first  gives,  on  substituting  for  v X,  X and  9, 

(a2q,  - qq,q)  * q,  = (2qq,  + <\t)3>  - + £ • q.,) 


(2.10) 


-X 
that  is 


-X 


.2  2 
a \ 


‘x  T na'-xt  ttvnx 

, o . „ 

+ r 


< “ (<wx  + \tr  - V*xx  + £ ' 2xx 


^XX' 

) (2.11) 


Since  the  ultimate  aim  is  to  replace  (x,y)  by  (X,9),  it  remains 
to  obtain  suitable  expressions  for  the  scalar  product  r * q 
occurring  in  (2,11),  This  is  achieved  by  noting  that,  except  in 
degenerate  cases,  the  original  field  equations  defining  X(x,y,t), 
®(x,y,t), 


$+r*q=9,  <J>^+r»q.=0 


(2.12) 


may  be  solved  explicitly  for  r to  yield  a parametric  representa- 
tion of  the  plane  in  the  form  r(X, 9,t).  The  position  vector  of 
the  field  point  occurring  In  (2.11)  may  be  replaced  as  required 
by  functions  of  the  independent  variables  (X, 9,t)  by  using  the 
expression 
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where 


q + A^(X)q  + B-i_(\)q  = 0 


Ax(\) 

B^M 


„3q  2 

= - q d 

2q  dX 
\ 


l* 


X * 


-4  -2 


- ^ log  q2ex 


(2.13) 


(2.14) 


Scalar  multiplication  of  (2.13)  by  r ^nd  using  (2,12)  gives 
finally 

r • q = (®  - (p)A1(\)  + B1(\) 

so  that  (2.11)  becomes 

a2£\  = (qq\  + ' ^tt^w  + \ B1  * A1(<T>  " ^ (2.15) 

o 

in  which  a , being  given  by  (2.10),  is  independent  of  9. 

<p(r;t)  never  reduces  to  a function  of  (X,t)  only  for 
otherwise  the  velocity  vector 


q = V?  = 9XVX 

will  only  give  irrotational  flow 


* = VA  = VXA  \ = ° 

in  the  trivial  case  q q = 0.  It  follows  that  if  the  field 

“A  \ 

equation  (2.15)  is  to  be  valid  everywhere  then 


A f\  4-  \ — r\ 

x to 

3.nd.  2 2 ^ 

a ax  = <mx  + + hV 

The  condition  (2.16)  shows  at  once  that  either, 
^^.(Xjt)  must  vanish  identically. 


(2.16) 

(2.17) 

or  both,  of  A-j_(X), 
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I*  (X,t)  = 0.  $j.(X,t)  reduces  to  a function  of  X only 

say  g(\)  and  we  can  write 

®(X,t)  = tg  ( X)  + f(X)  (2.17) 

where  g(X)  is  related  to  q(X)  by  (2,17) ? 

a2q2  = [qq^  + g‘(X)]S  (2.18) 

and  22  2 

%£■  + + 2^X)  - (2.19) 

It  follows  the  parameter  X(x,y,t)  is  given,  implicitly,  by  the 
field  equation  (2.4), 

tg‘(X)  + f»(X)  + £ • q = 0 (2.20) 

A 

Regarding  f and  g as  functions  of  q the  following  rela- 
tion is  equivalent  tc  (2.18),  and  (2.19), 

(Y  - 15(1  + q2eQ2)(c2  - q2  - 2g)  = 2(q  + g )2  (2.21) 

4 q 

The  functions  f(q);  g(q)  may  therefore  be  chosen  arbitrarily, 
the  velocity  distribution  0(q)  corresponding  to  any  selected 
g(q)  being  obtained  from  the  ordinary  differential  equation  (2.21). 
The  field  equation  (2.20)  may  then  be  used  to  determine  9(q)  as 
a function  of  (x,y,t) 


II.  A-^(X)  = 0.  The  second  possibility  implies  from  (2.13), 


and  (2.14) 


log  42«x)  = 0 


that  is 


d ( 1 _ 

^ ^ -X 


q J = o 
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so  that  (1/q  9^)q^  is  a constant  vector.  Therefore  the  X-lines 
(2.U-)  comprise  a system  of  lines  parallel  to  a fixed  direction. 
Integration  of  A^(X)  = 0 twice,  using  (k.l’1-).  gives 


q = q0  cosec  (©  + 9 ) 


(2.22) 


in  terms  of  constants  of  integration  q0,  9q.  The  component  of 
velocity  in  the  direction  9 =£“0O  has  the  constant  value  qQ, 
Since 

q sin  (9  + G)  = q0  t(9) 

with  the  usual  potation  for  t(9)  the  unit  tangent  vector  to  th£ 
streamline.  Then 

.2 


therefore  the  X-lines 


may  be  rewritten 


sin  (9  + 9q)  = - qQ  t(0Q) 


<p  + r * q = 
S S 


- O' 


2 ,• 


5>e  sin  (9  + 9q)  = q0  r • t(-  ©Q) 


(2.23) 


When  expressed  in  terms  of  9 the  equation  f2.17)  becomes 

= <qqe  + $et)2  - *tt  q2  H5 

and  then,  using  the  parameter  z = qn  cot  (0  + 9^),  may  be  trans- 

w V 

formed  into  the  form 

a2  = (z  + 0 )2  - $ $> 

zt  zz  tt 


wnere 


= 02  - q2  . . 2® 

Y~1  0 t 


or  finally  using 
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F(z;t)  = <J>  ( z , t ) + ^ tz 2 + ^ (q|  - c2)t  (2.2*0 

so  that  » 

a2  = - ( y - l)Ft(z$t)  (2.25) 

then  o 

F - F F.  . + tF  + . i)»  _ u (2.26) 

zt  m U tt  ‘ t 

Similarly  the  equation  (2.23)  for  the  parametric  lines  may  be 
expressed  in  terms  of  z and  F(z,t), 

F (z;t)  - tz  + x cos  9 - y sin  0=0 

z o o 

which  serves  to  provide  z as  a function  of  (x,y,t)wben  any  solution 
of  (2,26)  for  F(z$t)  is  inserted.  If  9 = 0 is  selected  then 

the  parameter  lines  are  parallel  to  the  y-axis  and  are  given  by 

x = tz  - F (z ; t ) 
z 

where  z = q cot  Q. 
o 

Concurrent  Parameter  Lines.  When  the  classical  theory 
of  the  steady  irrotational  simple  wave  is  developed,  the  general 
simple  wave  generated  by  X-lines  enveloping  a convex  curve  is 
obtained  by  a generalization  of  the  original  Prandtl-Meyer  expan- 
sion wave  for  which  the  X-lines  are  concurrent  through  a fixed 
point.  It  is  clear  however  that  this  procedure  cannot  be  repeated 
to  obtain  a generalized  unsteady  simple  wave  from  the  case  of 
concurrent  parameter  lines,  for  this  latter  is  the  very  case 
which  is  automatically  excluded  by  (2,17).  For  if  the  X-lines 
are  concurrent  through  a fixed  point,  this  may  be  selected  to 
coincide  with  the  origin  then  (2.20)  gives  g* (X)  = f * (X>  = 0 so 
that  0 reduces  to  a function  of  time  only,  in  fact  $ = at  t-  b 
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where  a,b  are  constants.  The  equation  for  the  velocity  (2»l8) 

reduces  to  ? 2 2 2 

a qx  = q q* 

where  2 2 

V + — a— — = constant 

Y - 1 

so  that  the  entire  flow  has  reduced  to  the  Prandtl-Meyer  simple 
centered  expansion  wave  in  steady  flow.  The  second  possibility 
A-^(X)  = 0 implying  parallel  X-lines  (preceding  pages),  it  follows 
that  there  is  no  unsteady  simple  centered  wave,  with  fixed  center, 
of  the  type  q(X), 


The  General  Case.  The  envelope  E(t)  will  be  obtained 
by  solving  for  r the  equations 


r 

+ L * = 0 

I 


(2.27) 


+ r • q..  =0  I 
XX  - — X\  ^ 

the  second,  representing  a one  parameter  family  of  straight  lines, 
being  obtained  from  the  first  by  partial  differentiation  v/ith 
respect  to  X (x,y,t  assumed  constant).  These  lines  will  be 
parallel  if  q q = 0,  or  A, (\)  = 0,  which  exceptional  case 
shall  be  excluded  at  present.  For  any  given  (X,t),  the 
system  (2.27)  represents  a line  pair  intersecting  in  the  point 
rQ(X;t)  given  by 


r (X;t)  = | q . q I 

o X-XX  XX  *X  1 — \a  -XX1 


-1 


Now  vector  multiplication  of  (2.13)  by  q gives,  since 

X 
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i 

q q J = q ©v  = exp  - 
“A  ~X 

then 

lq  q | s A (X)  exp 
— X> 

therefore 

A1(^)Ho(X; t)  = (<I>X  £xx  - $xx  J^)  exp  ] B1(X)dX  (2.28) 

o’ 

As  X varies,  rQ(X,t)  will,  for  fixed  t,  generate  a curve,  the 
envelope  E(t)  of  the  X-lines,  If  the  X-lines  are  concurrent  then 
E(t)  degenerates  into  a point  and  r^/X^t)  is  fixed  for  all  values 
of  X,  That  is  reduces  to  a function  of  time  only,  therefore 

9r  /3X  = 0. 

— o 

Elimination  of  the  vector  £ from  (2.13)  gives 

a\ 


B1(X)dX, 

| B1(X)dX 


Al(X)ro(X;t)  = - [$xx£x  + ^ (A^Q,  + B^)]  exp  j B1(X)dX 


giving  after  partial  differentiation  with  respect  to  X and 
repeated  use  of  (2,13), 


VM  -0\  ~ ' + Bl\>.  + (A1  + B1X)4>X 

A \ P 

- r-  (®xx  + B1  V]£x  x exp  I Bi(X)dX 

1 J 

so  that  (r  )\  = 0 when  OCXstl  satisfies 
~"0  \ 


J3_ 

ax 


(*XX  + BA)  - (Oxx  + B1$x) 


= 0 


giving,  on  partial  integration  in  terms  of  a function  of  inte- 
gration h(t), 
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Since  the  case  A^(X)  = 0 is  excluded,  then  <£(X;t)  is  restricted 
to  the  special  form  (2,17)  so  that 

(g\\  + Bxgx  + Axg)  + (f^  + B1fx  + A±nt  = A3(\)h(t) 

in  terms  of  g(X),  f" ( X) , A^(X),  B^(X)  which  are  not  functions  of 
the  independent  variable  t.  It  follows  that  h(t)  is  not  arbitrary 
but  must  be  a linear  function  of  t,  defined  in  terms  of  constants 
(aoJ  t>n),  not  necessarily  non-zero, 

h(t)  = a + tb_ 
o o 

so  that  fU),  g(X)  are  now  given  by 


EXX  + B1EX  + Vg  - ao>  = 0 


fxx  * Vx  + Ai<f  - V = 0 

Therefore,  apart  from  additive  constants  not  affecting  $x(X;t), 
the  functions  f(X),  g(X)  both  satisfy  the  ordinary  linear  differ- 
ential equation 


EXX  + BlgX  + Alg  ” 0 


(2.29) 


in  order  that  the  X-lines  may  be  instantaneously  concurrent. 

When  this  condition  is  satisfied  the  time  variation  of  r (t)  will 


— o 


be  given  by  (2,28);  since  = g(X), 


A.  ( X)  (r  ) = 

1 o t 


(Ex2xx  • Exx  V 9XP 


B1(X)dX 


or 


(r) 

o t 


“ \ + §xx  fixJ  e^p  | B1(X)dX 

B1(X)dX 


[g  fix  - gx  £]  exp 


• £0(t)  Being  independent  of  t or  X;  it  follows 


in  virtue  of  (2,29) 
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that  the  point  of  concurrency  (the  instantaneous  "center51  of 
the  wave)  moves  with  constant  velocity  in  a fixed  direction* 
[Although  apparently  a simple  function  of  X,  it  may  he  demon- 
strated by  direct  differentiation  with  respect  to  X and  the  use 

of  (2.29)|  (2.13)  that  (r  ) = 0,  as  it  must  since  (rQ).  = 0, 

0 X t A 

so  that  (rQ)  reduces  to  a constant,] 

Therefore  if 


$xx  + VX)<\  + vx)<&  = 0 


(2, 


«%  /*\  N 

ow 


the  X-lines  are  concurrent  through  the  point  rQ(t) 


Al(X)  £o(t)  = (Ox  0XX  - <^xA)  exp J ^OOdX 
which  with  the  aid  of  (2,13),  and  (2,30)  reduces  to 
rQ(t)  = (®fix  - ®x  fi)  exp 
which,  in  fact,  is  also  the  intersection  of  the  lines 


B1(X)dX 


$ + r • q = 0;  + r • qx  = 0 


(2.31) 


The  three  lines 


$ + r * q = 0;  4>  + r * q = O5  $ , + r • q.N  = 0 

are  therefore  concurrent,  for  all  X,  through  the  point  rQ(t). 

This  property  may  also  be  derived  by  a similar  method  to  that 
resorted  to  for  the  steady  rotational  flows;  equations  (2,13), 
and  (2,30)  being  exactly  analogous  to  (1,20)  and  (1,19), 

When  (2,30)  is  satisfied  the  unsteady  simple  wave  reduces 
to  a flow  which  is  essentially  "steady".  In  fact  a simple 
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transformation  to  the  moving  rectangular  coordinate  system  r 


i 


r‘(X,9,t)  = r(\,9,t)  - rQ(t) 


gives  since  r(X,9,t)  = [(&  - 9)^  - £.]  exp  J*  B1(X)dX 

obtained  by  solving  (2.3),  and  (2.4),  therefore 


r 1 ( X,  cp, t ) = - 9 0 exp 

\ 


B1(X)dX 


so  that  the  revised  representation 
equation  of  the  parameter  line 
moving  axes  is 

+ (r1  + r ) • 

A.  O 


reduces 

+ ^ • Sx 


to  r*  (X,9)  and.  the 
= 0 referred  to  the 


in  which  (2.31)  gives 

\ + ^ \ « • ad  exp  J 

= 0 

in  virtue  of  the  definition  (2,14)  of  B-^(X).  The  fundamental 
equation  of  the  lines  of  constant  X therefore  reduces  to 


r ' ♦ q = 0 

which,  since  q = q(X),  defines  X as  a function  of  r1  only,  not  t, 
Pef erred  to  the  new  coordinate  sys cem  moving  along  with  the  point 
of  concurrency  £ (t)  with  constant  velocity,  therefore,  the  flow 
reduces  to  a steady  phenomenon  - the  simple  central  Prandtl-Meyer 
expansion  wave. 
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Non-Degenerate  Envelope.  The  case  of  concurrent  X-lines 
representing  no  essentially  new  phenomena,  attention  will  be 
directed  to  flows  possessing  a non-degenerate  X-envelope,  For 
such  flows  the  functions  f , g may  be  regarded  as  functions  of 
q (this  amounts  to  selecting  the  arbitrary  parameter  X(xyt)  to 
coincide  with  q(xyt)  = q(X)  = q = X).  The  relevant  equations  are 

(Y-  D(1  + q^©q2)(c2  - q2  - 2g)  = 2(q  + gq)2  (2,32) 

tg  + f +r*q  =0  (2.33) 

q q -q 

then,  provided  at  least  one  of  the  relations 


Snn  + + A1 Cq)g  * 0 

qq  -l  q 1 

f + B (q)f  + A (q)f  * 0 
x q j- 

is  not  violated,  th^  X -lines  will  not  be  concurrent, 

2 

Solution  1.  Take  for  example  g(q)  = i-  q f(q)  = 0 then, 

if  b2  = c2/2  0 

a2  = ( Y ~ 1) (b2  - q" ) 

let  p p 2 2 ? 

q = b (sin  z + p cos  z) 


then 


p2  = Ur. _i 

Y + 3 
a = 2bp  cos  z 
so  that,  from  (2,32), 

o -1  x 

qQ  = p tan  z 

q K 


(2.3*+) 


when  the  positive  sign  for  the  square  root  is  selected. 


Expressed  in  terms  of  z, 
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p@  = (1  - p^)  &+p^cot2z)  ^ 
z 

which  may  be  integrated  to  give 

© = p ^z  + tan  ^(p  cot  z)  + n - — (2,35) 

2p 

by  suitable  choice  of  the  constant  of  integration. 

Now  from  (2.5)  the  angle  which  the  streamlines  make  with 
the  lines  of  constant  velocity  is  given  by. .using  (2,34), 

tan  P = p cot  z (2,36) 

so  that  (2.35)  may  be  rewritten, 

© = p"1*  + P + (2  - p"1)  l-  (2,37) 

(2.33)  becomes 

qt  + x(q  cos  ©)  + y(q  sin  ©)  =0 

or,  on  performing  the  differentiation,  using  q©^  - cot  P;  and 
expressing  x,y  in  terms  of  polar  coordinates  (r,cp)  (x  = r cos  cp  , 
y = r sin  9) , 

qt  sin  P = r sin  (S-P-cp) 

Finally,  in  terms  of  z,  using  (2,36),  and  (2.37), 

R sin  (p-1z  - 9 - -£•)  = cos  z (2,38) 

2p 

where  R = r/pbt.  For  any  fixed  z,  equation  (2.38)  represents  the 
(R,9 )equation  of  a straight  line  and  the  envelope  of  the  system 
as  z varies  will  be  given  additionally  by  partial  differentiation 
with  respect  to  z;  yielding  a second  line, 


f 


i 


* 


I 


\\ 


j 


jr\ 

i 

m 

e 

B 


All -111 


"^1  TC 

R cos  (p  z - cp  - £~)  = - p sin  z 


26 

(2.39) 


and  therefore  the  equations  in  parametric  form  may  be  expressed 
in  terms  of  z (in  terms  of  which  T,  q,  6 are  already  known). 

<p  = p ^z  - tan-'1'  (p  tan  z)  - ~ + ^2! 


2p  T 


2 2 2 9 

R = cos  z + p sin^z 


(2.40) 

(2.41) 


therefore  the  (R,cp) 
desired,  since 


equation  may  be  explicitly  written  down  if 


tan  z = (_L^X)1/2 


•'2  .2 ' 
r - p 


and  0 < p < R < 1 (p2  = -L  for  Y = J) 

- - ip  '5' 

If  In  (2 , 40) , z is  regarded  as  the  radial  distance  in 
polar  coordinates  (z,cp)  then  (2,40)  is  the  equation  of  an  epi- 
cycloid, and  since  (2,4l)  may  be  rewritten 


R2  •-  (1  - p2)  cos2z  + p2 


the  actual  curve  in  the  (R,<p)  may  be  drawn. 

The  envelope  is  similar  to  an  epicycloid  which  increases 
linearly  in  size  since  r = R pbt.  Along  each  tangent  line  the 
velocity  vector  is  instantaneously  constant  and  is  inclined  to 
the  Dine  at  the  angle  P according  to  (2.36),  The  instantaneous 
streamlines  are  perpendicular  to  the  tangent  generated  from  the 
point  of  the  envelope  given  by  z = 0 and  then  sweep  round  to 
approach  asymptotically  the  direction  parallel  to  the  tangent  at 
z = it/ 2 • 

The  parameter  line  z = 0 along  which  the  flow  commences 


r 


■a-  • 3 
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x 


moves  normal  to  itself  with  velocity  pb,  equal  to  that  of  the 
fluid  along  it.  The  z = 0 line  is  perpendicular  to  the  radius 
© = - P where,  p - - n/2 (p~^-3 ) lies  in  the  fourth  quadrant, 

(0  < p < ti/2  for  y = 7/5;  p2  = 1/11). 

3 • Solution  Involving  an  Initial  Advancing  Discontinuity. 
When  the  X-lines  are  parallel  the  flow  is  given  by  a solution 
F(z;t)  of  the  non-linear  partial  differential  equation 

Flt  - F22  Ptt  + * Ptt  + (Y-  1)Ft  0 (3a) 

the  velocity  parameter  z(x,t)  being  then  given  by  the  field 
equation 

x = zt  - Fz(z;t)  (3.2) 

For  the  acoustic  speed;  a “ = X_rJL_(c2  - q^  - - 2$’^) 

that  is  a^  = - (y  - l)Ft(z;t) 

so  that  for  non-negative  pressures  we  must  consider  solutions 
of  (3.1)  for  which  the  time  derivative  is  negative  only. 

An  interesting  solution  exhibiting  an  unsteady  singular- 
ity in  two-dimensional  flow  may  be  obtained  by  separating  the 
variables  and  trying 

F(z ; t)  = Z(z)T( t ) 
giving  as  a possible  solution, 

tTtt  + - nit  = o 

provided 

7Z  t2 

— = ~ — = constant,  independent  of  z or  t, 

T_tt 
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The  first  implies  T(t)  = A^t2^  (in  terms  of  a constant  A)  and 
this  is  automatically  compatible  with  the  second  and  third 


conditions  when 


Z(z)  = A2z' 


so  that  we  may  take  as  a special  solution  of  (3,1), 

F(z’ t)  = - ^ Y-l  2-Y 

F(  J } (Y  - 1)(2  - Y)  t 


(3.3) 


by  suitable  choice  of  the  constant  b which  must  be  positive  for 
1 < y < 2.  Taking  y = 7/5  then  the  parameter  z(x,t)  is  given 


implicitly  by 


« Zt  + (i)3/5 

3 z 


(3.4) 


f z > 0 i 

Confining  attention  to  ^ ^ ">  0 J ’ this  aquation  defines  z real 


only  when  x exceeds  its  value  when 


- zxU ) - 


at  which 


= ° 
5z 


which  implies  a concentration  of  differing  z lines  (wave  lines) 
and  rapid  change  in  velocity  at  the  station 

X = | b5/8t3A 
= x1(t)  say. 

The  wave  line  specified  by  constant  z advances  with  speed 


If  = Z 4-  bz-3/5t-2/5 

which  is  very  large  near  t = 0 and  decreases  to  a constant  -value 
z,  equal  to  the  x-component  of  the  instantaneously  local  fluid 
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speed,  as  t->oo  , In  contrast,  the  singular  line,  which  at  t = 0 
is  crossing  the  y-axis  with  00  velocity,  decelerates  with  speed 
eventually  falling  to  zero  as  t — ^00  . This  speed  is 

A (|  b?/8t3/S  = 2b^/8t_1/lf  = 2z1(t) 

is  twice  the  x-componer.  t of  the  local  fluid  located  on  the  line 
at  time  t.  The  instantaneous  streamlines  are  defined  by  dy/dx  = 
qQ/z  so  that,  regarding  t as  constant, 


bt  , 

dz  z dz  " z ^ £75' 

z 


therefore 


(y  - y0)/qc  = t log  z + t.3//5z-8/5 

and 

x = tz  + 2*  t3/:?z~3/5 

3 

may  be  regarded  as  the  parametric  equations  of  the  streamline 
in  terms  of  z as  parameter,  yQ  being  a constant  determining  the 
different  streamlines.  Singular  points  in  all  such  streamlines 
will  occur  when  3y/3z  ~ 9x/3z  = O5  but  the  curvature  is  infinite 
and  dy/dx  remains  finite  and  one  valued.  This  occurs  when 
z = z^(t)  and  it  is  seen  that  no  streamline  may  be  continued 
beyond  the  singular  line  into  the  region  x_<x^(t),  corresponding 
to  which  no  real  values  of  z(x,t)  exist. 

Also  both  branches  of  the  x~z  curve,  corresponding  to 


the  variations 


0 < z < z-^(t) 


zx(t)  < z <00, 


t+f  . 
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cannot  be  simultaneously  considered  at  a given  positive  time  t, 
because  this  would  imply  intersections  of  the  streamlines  and  a 
double  valued  velocity  at  the  point  (x,y,t). 

However,  in  unsteady  motion  the  instantaneous  stream- 
lines do  not  yield  a picture  of  the  actual  flow  pattern  and  it 
is  necessary  to  obtain  explicitly  the  equations  of  the  paths  of 
the  fluid  elements.  These  are  obtained,  theoretically,  by 
integrating  the  equations, 


dx  _ z/x  dy  _ 

at  3T  qo 

the  parameter  z being  given  by  (3.4).  It  is  necessary  to  adopt 
an  alternative  procedure  and  to  eventually  express  x and  y in 
terms  of  z,  not  the  time.  Form  (3*4)  by  taking  the  total  differ- 


ential 


dx  = |$  dt  + dz 
Ot  Oz 


therefore  the  relation  between  z and  t along  a particle  path  is 
given  by  the  ordinary  differential  equation,  for  z(t); 

3x  , dx  dz 

at  Sz  dt 


that  is,  using  (3.4)  explicitly; 

. .7/5  3/5 

dt  _ t t z 
dz  z b 

giving  an  integration  in  terms  of  a constant  a,  depending  on  the 


element, 


p 2/5 

2 *■  ao  = 5b(f) 


(3.5) 


Therefore  on  substitution  for  t as  a function  of  z; 
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4z2  + a 
(z2  + a)?/2 


y - yr 


q0(5b)- 


(z2  + a)' 


(3.6) 


are  the  equations  of  the  particle  paths  in  which  (ajqo)  determine 
2 

the  oo  elements  and  z,  the  x component  of  the  fluid  speed  is 
used  as  fundamental  parameter.  The  general  form  of  the  fluid 
element  paths,  according  to  (3.6),  will  differ  according  as  a 
is  positive  or  negative. 

It  is  seen  that  when  a > 0 the  fluid  paths  possess  a 
cusp  at  z = 1/2  • Then  the  element  (a,qo)  will  generate  this 

cusp  at  time  t = 16  b^2/a2  along  the  line  x = 64/3  b'^cc  ~>^2 
that  is  when  it  meets  the  line  x = 8/3  b^^t^^  i.e.  x = x^(t). 
When  this  occurs,  the  local  acoustic  speed,  is  given  by 


a2  - b Cf )2/?  = l 


('1.7') 


so  that  I a I = ^ = z-^(t).  Therefore,  on  the  singular  line, 

the  streamlines  intersect  this  line  at  the  local  Mach  angle 
“ sin”^  |S|  ; or  the  normal  component  of  the  fluid  speed  is 
sonic.  In  particular  it  follows  that  the  singular  line  only 
occurs  when  the  local  fluid  speed  (q2  + z2)  > a^  is  supersonic, 


The  acceleration 


7 4-  n • VJ  *7  — *7  4-  */  7 

zt  ■ 2.  v “ “t  x 


= bz(z8/^V/?  - bt)”1 


5/8  -1/4 

becomes  infinite  when  z = b t = z^(t)  so  that,  in  effect, 
the  solution  breaks  down  on  and.  to  the  left  of  the  line  x=x^(t), 
Any  attempt  to  continue  the  flow  into  the  region  0 < x < x1(t), 
which  was  previously  occupied  by  fluid  in  continuous  motion, 
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must  fall  because  the  field  Eqn.  (3.4)  will  not  give  a real 
velocity  for  x < x^(t).  The  parametric  equations  (3*6)  are 
similarly  not  applicable,  and  cannot  yield  real  z for  real 
x < x.,  ( t ) « 

X 

First  Solution,  When  elements  a > 0 are  considered  there 
are  two  possible  patns.  The  first  commencing  at  t = 0; 

c* /p  — ”2  /p 

x = 4/3  (5h)'v  a-'5  continues  to  move  to  the  right  until  it  is 

5/2  -2 

eventually  overtaken  at  time  t = 16  b a by  the  singular  line 
x = x-^(t)  which  by  that  time  has  reached  the  station 
x = 64/3  b^/2a“-^//2[  = x1(t)].  The  relevant  range  for  z is 
0 < z < 1/2  s/a  . When  a < 0,  the  corresponding  fluid  elements 
are  concentrated  along  the  axis  x = 0 at  zero  time.  This  diffi- 
culty avoided  by  confining  attention  to  a > 0 for  which  at 
t = 0 the  instantaneous  streamlines  occupy  the  whole  half  plane 
x >_  0 and  comprise  lines  parallel  to  the  x = 0 axis  since  z = 0, 
The  singular  line  initially  found  along  x = 0,  moving  at  high 
speed,  may  be  regarded  as  the  result  of  a uniform  impulsive 
pressure  applied  over  this  face  (the  x = 0 plane  in  three-dimen- 
sional flov ) from  which  compression  waves  are  generatedj  These 
waves  proceed  outwards  with  ever  decreasing  velocity,  deflecting 
the  fluid  particles  from  their  initially  vertical  paths. 

The  compression  waves  precede  the  advancing  singular  line 
x = X-|  (t ) which  overtakes  the  fluid  instantaneously  moving  where 

u. 

x-component  of  velocity  is  locally  sonic.  At  the  instant  of 
interception  this  component  is  abruptly  doubled  to  equal  the 
velocity  2 z^(t)  of  propagation  of  the  discontinuity,  the  fluid 
being  carried  along  in  an  ever  increasing  concentration  (density) 
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along  the  singular  line.  This  process  causes  the  speed  of  pro- 
pagation to  reduce  continuously  and  fall  to  zero  as  time  increases. 
It  should  be  noted  that  the  advancing  singular  line  cannot  be 
considered,  as  a shock  wave;  the  fluid  overtaken  has  its  speed 
discontinously  increased  but  the  normal  velocity  relative  to 
that  of  the  front  rereduces  to  zero  so  that  this  change  does  not 
resemble  passage  through  a shock.  In  addition  the  fluid  does 
not  pass  through  the  singularity  but  accumulates  along  the  front 
and  is  carried  along  with  it,  leaving  a vacuum  in  the  region 
behind, 

7 ? 2/5 

Since  p £a'  and  a = b(~0  , the  uniformly  distributed 

surface  pressure  necessary  to  be  applied  along  the  boundary 

5 78  -l  /b 

7.  - 2C,  (t)  is  obtained  when  z = zn  (t)  = b t therefore 

_ r +-iA 

P Q_  t o 

Second  Solution.  When  the  second  branch  z^(t)  < z < go 
is  considered  flow  obtained  is  quite  different  because  at  t = 0 
there  is  an  intense  concentration  (hjgh  density)  of  fluid  along 
the  x = 0 axis,  across  which  it  is  instantaneously  released  with 
high  speed,  expanding  into  the  adjacent  vacuum  x > 0 (x  < 0 can 
be  similarly  considered).  At  a subsequent  time  t the  singular 
line  has  advanced  to  x = x-^(t)  with  ever  decreasing  speed  and  the 
fluid  released  at  t = 0 continues  to  precede  the  singular  line 
also  with  ever  decreasing  speed.  Unlike  the  phenomenon  for 
a > 0,  the  front  never  overtakes  the  preceding  a<0  fluid  again 
but  falls  back  as  time  progesses.  This  may  be  demonstrated  by 
focusing  attention  on  conditions  at  a chosen  time  t(>  0).  At 
this  time  the  front  is  moving  forward  with  speed 
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v = 2z-^(t)  = or  t = 2lfb^//2v  \ and  the  x-component 

of  velocity  of  any  particle  a < 0 is  z given  by  t=(5b)^^z(z2-Ki)  2 
5/2  „h. 

_>  (?b)  z • Therefore 

>+  5/?  _k  5/2  -*+ 

?.  b5/"v  > (5by  z 


.5/8-1 


z > 5 2 V > V. 


It  follows,  since  this  inequality  is  valid  for  all  positive  time, 
that  all  the  a < 0 fluid  elements  initially  concentrated  there 
precede  the  singular  line.  A similar  argument  shows  that  the 
a >0  elements  are  characterized  by  z < v,  so  that  considerable 
information  is  obtained  from  the  curve  showing  the  variations 
of  x i^ith  a (-oo  < a ^hb^^^t-1^2 ) at  any  tim~  t > 0. 

This  may  be  obtained  by  noting  that  >:(a,t)  is  given 
jointly  In  terms  of  (z.a)  by 


4z2  + a 
(z2  + a)5/ 


(3.8) 


where  z(a,t)  is  given  by 


(z2  + a 


(3.9) 


Ir  follows  by  differentiation  partially  with  respect  to  a 
(t  constant)  for  x(a,t)5 


3Z  g,t(4z2-0)|f  ♦ * 

2 ( H-z2  - a)  p.  + 5z  = 0 
oa 


Therefore,  eliminating  Sz/9c., 
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0X  t_ 

0a  2z 


(3.10) 


The  subsidiary  curves  giving  the  variation  of  x,a  with  z are 
first  drawn  and  then  combined  to  give,  with  the  aid  of  (3.10), 
the  curve  x.  ^ a • These  curves  have  equations 


x = zt  + 


2k  {t) 

3 V 


8/5 


a 


= 5b  <§> 


2/5 


- z 


and  generate  a cusp  at  [x-i(t),  a-^(t)]  in  the  (x,a)  locus  where 


x 


(t)  = | b?/8t3A;  a1(t)  = At“1//2 

_ K5/8+-lA 


1' w/  _ 3 “ w » ~l 
„(t)  = (5b)5/8t"lA; 


5 zx(t)  = 


The  curve  varies  slightly  with  t but  it  is  seen  how  only 
one  branch  may  be  considered  at  a time,  since  for  any  given  (x,t) 
two  different  values  of  a , with  corresponding  different  velocities 
cannot  be  simultaneously  considered.  As  t increases  the  cusp 

v ^/li  /li- 

CX  = H-h"  ’t  * approaches  zero,  the  total  variation  of  cx  xn  the 
first  solution  [0  < a _<  M-b^^f-^2 ; 0<  z<_  z-^t)]  being  con- 
tinuously reduced  as  more  and  more  elements  are  caught  in  the 
advancing  front.  A similar  state  of  affairs  exists  when  the 

if  /\x  -1 

second  solution  [ -co  < a <bo'  t~  zi_(t)  < z < + oo]  is  con- 
sidered, particles  for  which  a > 0 being  removed  from  the  regime 
of  continuous  fluid  flow  by  the  ad  Tancing  front,  but  the  a _<  0 
elements  are  never  reached  and  continue  their  motion  uninterrupted 
for  ail  time. 


The  variation  of  speed  with  time  specified  by  (3.9)  of 
elements  for  any  a ; -oo  < a < + oo  relative  to  that  of  the 
singular  line  which  is  z = 2z]_(t)  = 2b^8t~**’//1+  shows  that 
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1.  Elements  a < 0 always  .move  faster  than  the  singular  line 
at  any  time  t. 

2,  Elements  a > 0 in  range  \/a  <_  z < oo  initially  move 
faster  but  are  eventually  overtaken  by  the  singular  line 

These  results  follow  because  the  intersection  of  the 
singular  line  curve  A with  the  a > 0 curve  B occurs  in  the  range 


1 

2 


< z < I a 


1/2 


The  front  of  the  disturbance  which  expands  into  the  vacuum  is 
therefore  composed  of  high  speed  a < 0 elements. 

The  field  equation  curve  is  shown  in  Fig,  1.  The  Figures 
3,  4-,  and  6 are  used  to  draw  the  particle  paths  and  Fig.  2 the 
instantaneous  streamlines  of  the  first  solution.  The  correspond- 
ing curves  for  the  particle  paths  of  the  second  solution  are 
shown  in  Figures  7,  8,  9,  and  10. 


4,  Generalization  of  the  Unsteady  Simple  Wave.  When 
the  direct  extension  of  the  simple  wave  q (\)  was  considered  in 

Section  2,  the  combined  single  equation  for  the  velocity  (obtained 

\ 

p 

by  simple  elimination  of  a . between  (2.1)  and  (2.2))  yielded  the 
two  equations  (2.16),  (2,17)  which  could  be  replaced  by  an 
ordinary  differential  equation  for  the  velocity  together  with  a 
partial  differential  equation  for  0(\,t).  There  were  two  dis- 
tinct possibilities  but  both  were  restrictive.  Either  the  \-lines 
were  parallel  (in  which  case  the  flows  could  be  treated  by  a 
more  direct  method),  or,  in  the  case  of  a non-degenerate  envelope, 
the  time  variation  of  3>(\,t)  was  linear.  If  the  preceding  method 
is  adapted  to  flows  of  the  type  q(\,t),  in  which  the  time 
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variation  occurs  independently  as  well  as  through  the  parameter 
X(x,y,t),  by  replacing  (x,y,t)  by  (X,9,t)  as  fundamental  indepen- 
dent variables,  then  the  original  single  partial  differential 
equation  for  the  velocity  regarded  as  a function  of  the  three 
independent  variables  (x,y,t)  may  be  reduced  to  three  partial 
differential  equations  for  q,  9,  regarded  as  functions  of  the 
two  variables  (X,t).  As  in  the  simpler  one  parameter  velocity 
field,  this  reduction  is  achieved  by  repeating  the  two  fold 
transformation  of  (x,y,t)  into  (X,9,t),  simultaneously  using  the 
Legendre  contact  transformation  to  transform  the  dependent  func- 
tion 9 into  O,  whilst  retaining  9 as  a fundamental  independent 
variable. 


The  reduction  achieved  is  again  due  to  the  fact  that 
*(x,y,t)  reduces  to  a function  of  (X,t)  only  and  is  independent 
of  9.  Such  flows  for  which  the  velocity  vector,  reducing  to  a 
function  of  two  independent  variables  (X,t)  are  sometimes  referred 
to  as  "double  waves"  but  the  title  simple  wave  is  retained  for 
the  flows  here  discussed  because  the  hodograph  in  the  velocity 
plane  still  reduces  to  a curve,  though  this  curve  may  vary  vtith 
the  time  and  generate  a surface  in  the  velocity  time  space® 

The  original  analysis  may  be  repeated  by  remembering  that 
since  q = q(\,t), 


aa 

■St  " a 


ax  . n 
x at 


in  which  q does  nor  mean  the  total  time  variation  of  q(X,t)  but 

“t 

simply  the  partial  derivative  with  respect  to  time  when  X is 
constant.  In  the  following,  therefore,  the  suffix  t is  used  to 
denote  partial  differentiation  of  a fiinction  of  the  type  f(X,t) 
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with  respect  tc  time  when  X is  constant  and  not  used  to  refer  to 
the  total  differentiation  with  respect  to  time  for  which  (x,y)are 
constant,  except  for  9 and  X which  are  regarded  as  functions  of 
(x,y,t).  The  following  revised  equations  are  obtained  by  successive 
differentiation 


0 + r 

• a = 

T 

(4.1) 

\ * E 

* % = 

0 

(4.2) 

ft  = °t 

+ r . , 

C^t 

(4.3) 

*tt = 

®tt  + e • att 

+ (\t 

+ L 

• q )X 
*Xt  t 

(4.4) 

(\x+  E • 2XX>*X  + . 

ix  = 

0 

- 

((J> 

XX 

+ E - ^xx)xt  + 

t4,xt  + 

r • 

a>t}  = 0 

whilst  the  equations  (2.10)  for  the  velocity  must  be  replaced  by 


(a2q^  - qq^  q)  * 7X  = 2q(q^Xt  + qt)  + $ 


V + _£ + *•  r • q — 

Y - 1 t - 2 

which  becomes  after  substitution  for  and  VX,  X , 

uX  J X/  * 


(4.5) 


(a2q  - qq  q)  * q = (2qq  + + r • a 1 ) (<D  + r • a.  ) 

—X  A — “X  A.  At  —\t  \t  — Xt 

“ (2qqt  + *tt  + £ ’ + £ 


(4.6) 


It  is  noted  that  9^,  9 , an<^  must  continue  to  be  regarded  as 

general  functions  of  position  and  time  and  the  object  is  now  to 
find  their  equivalent  representations  as  functions  of  (X,cp,t). 
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(4.7) 


The  required  expressions  analogous  to  (2.13)  are 

i \ 

+ + = 0 

3-Xt  + A2(X}t'I  + = 0 

q + A..  (X,t)q  + B.(X,t)q  = 0 

“tt  J - i -X 

I 

q + A,  (X,t)q  + B4(A,t)q  = 0 j 

t A J 

where  the  coefficients  A_^,  B_^  are  functions  of  (X,t)  only  and 

are  given  by  relations  of  the  type 


A,  q q = q q 
1 “A“X  “XA  “XX 


B1  q q = - q q 
1 “A~X  “A“XX 


therefore  explicitly 


"1 

B,  = 


A2  = 


3 2 
* °X 

a 

2qx 

ax' 

- wc  a°s 

& 

qx  + q2®x 

a 

ax 

qx  + q2^x 

a 

at 

q©£  0tt  + 

< 

" qte>. 

^X 

-1  q®x 


“2 » 

q 9X 


qe, 


(4.8) 

qe. 


B„  = 


B = - 


-1 

q2g  at 
q ^x 


e. 


x 
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The  actual  computation  of  these  coefficients  is  best  done  as 
follows.  With  the  usual  notation  for  t_(©)  = (cos  ©,  sin  S) 
the  unit  tangent  vector  along  the  streamline, 

q = qt (©) 

therefore 

= Q t + q©  n 

A.  ?S  ^ 

%x  = Cqxx  ■ ^x^  + (2qxex  + q9xx,a 

therefore 


or 


where 


q\qXX  + PCqxexx  - ex<‘u)  + 2qxe  + q2#xJq 

= (2qx9x  + q9xx)q% 

q^  + A-j^Xjt)^  + B1(X,t)q^  = 0 

q^©&l(X,t)  = (q2  + q2©^2)©^  + q2(— — )^ 


A 


after  some  rearrangement, 

2_  ......  ..2  2.2  N S 


^©^  A-^(X,t)  = (q^  + q^9^)  ~ E©  + tan-1  ~ ~ ] 


or, equivalently, 


r.  /%  \ 


a . . -2  <£ 

2qx 


Similar  procedures  yield  the  explicit  expressions  for  the 
remaining  coefficients. 
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Scalar  multiplication  of  equations  (4.7)  by  r and  then 
using  (4.1)  and  (4.2)  to  substitute  for  the  scalar  products 
1 • q 5 r • q,  gives 


r • qxx  = (4,  - Cp)A1  + 4^ 
r • q - (<S>  - cp)A2  + OxB2 


l 


r * q = (O  - <p)A-,  + <I>.B 
“tt  -3  ^ j 


(4.9) 


r • q - (0  - cp)Al+  + <&  | 

and  direct  substitution  into  (4.6)  and  (4,5)  gives 


X(\,t)  + (4>  - cp)Y(X,t)  + (4>  - q>)  Z(X,t)  = 0 


(4.10) 


in  which  the  functions  X,Y,Z  are  obtained  as  functions  of  c[(X,t), 
$(X,t)  only  and  are  therefore  independent  of  cp . It  follows  that 
if  the  field  equation  (4.10)  is  to  be  valid  everywhere  then  all 
the  coefficients  must  be  illusory,  tnat  is 


X(X,t)  = Y(X,t ) = Z(X,t)  ~ 0 


(4.10a) 


which  are  the  three  required  partial  differential  equations  for 
q,  ©,  O in  terms  of  the  independent  variables  (X,t), 

Explicitly 


Z(X,t)  = a2  - axa3 


(4.11) 


which  is  independent  of  <5(X,t) 


All-111 


42 


X(X,t)  = (®xt  + 2qqxx+  + B^) 

- <4>tt  + 2«t  - h\H\S  *l\>  - (Y  ' 

- \ - \\  j ix  + q2q^  (4,12) 

Y(X,t)  = 2A2(*xt  + qqx  + B^)  - ^ C ^ + B^) 

- A1  (<5tt  + 2qqt  + B^)  + ( y - D\  q?  (4.13) 


In  general,  solutions  of  the  equations  (4.10a).  (4.11),  (4.12), 

and  (4.13)  are  very  difficult  to  obtain,  such  flows  corresponding 

to  the  general  double  wave  solution.  The  original  simplified 

theory  of  the  unsteady  simple  wave  may  be  extracted  by  noting 

that  in  that  case  all  the  derivatives  q+ , ^ , q vanish  there- 

\t  —tt 

n "R  “ T3  — _ D c r>  4-  V»  vi  • 

0 4-  2 3 q 


fore  A2  = A_,  = A^  = B0  = B.,  = B,.  - 0 so  that 


X(X,t)  = (qqx+  4>tt)2  - 4>u(  \x  * B^)  - a2q2 

Y(X,t)  = - A,  ( X)4> 
x tt 

Z(.x,t)  = 0 


The  condition  Z = 0 is  therefore  automatically  satisfied  and  the 
additional  relations  X = Y = 0 reproduce  the  equations  (2,16), 
and  (2.17). 


The  X-Envelone.  When  the  generalization  of  the  simple 
wave  to  steady  rotational  compressible  flow  it  is  found  that  the 
X-lines  are  concurrent  or  parallel  (which  may  be  regarded  as  a 
particular  case  of  concurrent  lines  when  the  point  of  concurrency 
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recedes  to  infinity)  and  that  in  no  case  is  a non-degenerate 
envelope  possible.  This  is  in  contrast  to  the  original  steady 
irrotational  simple  wave,  for  which  an  arbitrary  curve  may  be 
assigned  in  advance  to  be  che  envelope  of  the  X-lines  from  which 
the  flow  is  constructed.  In  rotational  and  unsteady  incompress- 
ible flow  it  is  also  found  that  there  is  no  non-degenerate 
envelope.  The  simple  wave  in  "unsteady  flow  however  represents 
a much  more  general  class  of  flows  because  it  is  possible  to 
ootain  solutions  for  which  the  X -lines  are  neither  parallel  nor 
concurrent.  The  envelope  E(t)  would  be  obtained,  theoretically, 
by  solving  for  r the  equations 

1 


O.  + r 
K 

\\ + r 


^ - ° j 


(4.14) 


This  would  be  possible  provided  cu  q ? 0 (i.e.  A ^ 0),  the 

XA— XX  x 

exceptional  case  occurring  when  the  X-lines  are  all  parallel. 


If  the  lines  are  concurrent , then  the  point  of  concurrency  may 
vary  with  the  time  as  would  the  gens  al  envelope  vary  with  the 
time.  If  however  they  are  concurrent  through  a fixed  point, 
this  pejnt  may  be  selected  to  coincide  with  the  origin  without 
loss  of  generality  and  then  <5.  (X,t)  - 0 so  that  * reduces  to  a 
function  of  time  only. 

To  construct  the  flow  it  is  necessary  to  draw  the  envelope 
which  is  given  in  parametric  form  r(X,t)  from  equations  (4.14). 

On  the  envelope  there  is  a relation  between  X and  9 given  by 


\x  + ($  " CP)A1  + \B1  = 


0 


(4.15) 
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in  virtue  of  (4.14)  and  (4,9),  it  follows  that  v'v  and  X^  become 
infinite  on  E(t)  and  that  the  acceleration 


q (X,y,t)  + (q  • V)a  = (Xt  + q * ^ 

also  becomes  infinite.  The  solution  cannot  be  continued  as  far 
as  the  envelope  where  the  X-lines  coalesce  and  the  velocity  is 
ambiguous.  However  the  flow  field  of  any  solution  is  obtained 
automatically  by  drawing  the  X-lines  generated  as  the  tangents 
to  the  envelope.  Along  any  such  X-line,  the  value  of  X and 
therefore  of  q(X,t)  is  constant  and  known  from  (4,l4)» 

The  parametric  representation  of  the  plane  is  obtained 
from  (4.1),  and  (4.2), 


r (X,cp,t ) = [(<&  - cp)£x  - exp 


B1(X,t)dX  (4.16) 


since  q^9^  = exp  -J  B^(X,t)dX,  and  £ = By  direct  differ- 

entiation 


r = - + W*  A1(<|>  ~ cp)  ] 2 exp  f B^(X,t)dX 


r^  = - expj 


• v - \ 
O.A 


The  relation  between  X and  cp  along  an  instantaneous  streamline 
will  be  obtained  by  integrating  the  differential  definition 


q .dr  ~ 0 where 


dr  = r dX  + r dep 
“A  cp 


therefore 


+ A^(C>  - cp)]qdX  + q^  dep  = 0 


In  unsteady  flow  however  it  is  the  actual  paths  of  the  fluid 
elements  which  show  the  motion.  Those  paths  are  obtained  by 
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integrating  the  differential  equations 

r(X,<p,t)  = q(X,t) 

that  is 

ft  + *,  ft  ♦ *t  * *<»••«  (4'17) 

from  which  are  obtained  two  ordinary  simultaneous  differential 
equations  for  X and  9 regarded  as  functions  of  t only. 

On  integration  in  terms  of  two  constants  (c-,,Cp),  which 
serve  tc  determine  the  initial,  position  of  the  cc^  elements,  X 
and  9 are  obtained  as  functions  of  time  which  when  substituted 
in  (4.16)  yield  the  position  vector  of  the;  element  (c-pC2)  in 
the  form  rCtjC-pCg)* 

Singularities.  An  inherent  disadvantage  of  inverse 
methods  such  as  the  one  here  employed  is  that  iu  is  difficult  to 
avoid  in  advance  solutions  of  the  final  equations  giving  rise  to 
discontinuities  in  the  streamlines  in  the  physical  plane.  The 
velocity  may  be  a many  valued  function  of  position  corresponding 
to  those  regions  in  the  physical  plane  for  which  the  solution 
may  not  be  mathematically  unique,  but  it  is  possible  to  obtain 
a physically  possible  flow  pattern  by  confining  attention  to 
just  one  branch  of  the  X ^ r curves.  In  the  physical  plane  the 
instantaneous  streamlines  and  particle  paths  corresponding  to  the 
different  branches  usually  meet  in  a cusp,  where  the  acceleration 
becomes  infinite  and  the  flow  is  supersonic.  Each  solution  can- 
not bo  continued  beyond  the  cusp  locus  (limit  lino)  or  even  up 
to  it  because  such  singularities  represent  a local  breakdown  of 
the  continuous  flow  tneory  on  which  the  flow  was  constructed. 
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The  fact  that,  the  acceleration  becomes  infinite  where  the  local 
fluid  speed  is  supersonic  suggests  that  it  should  be  possible  to 
replace  the  limit  line  by  a shock  curve,  however  it  is  not  usually 
possible  to  satisfy  all  the  necessary  shock  relations  along  the 
limit  line.  In  the  simple  one  parameter  field  q(X)  of  Sections 
2,  3,  this  substitution  was  precluded  because  the  limit  lines 
intersected  the  streamlines  at  the  local  Mach  angle,  and  on 
existing  shock  theory  any  shock  inclined  at  the  local  Mach  angle 
has  vanishing  intensity  simply  reducing  to  a Mach  wave. 


The  Simple  Centred  Wave.  The  equations  (4-.  11),  (4. 12), 
and  (4,13)  are  so  complex  that  it  is  desirable  first  to  investi- 
gate the  possible  existence  of  flows  for  which  the  X-lines  are 
concurrent  through  a fixed  point,  and  defer  the  general  study  of 
these  equations  until  Section  5.  There  is  also  the  additional 
possibility  that  the  point  of  concurrency  may  vary  with  the  time. 
In  the  simpler  case,  however,  the  origin  may  be  selected  as  the 
fixed  point  through  which  all  the  X-lines  pass  for  all  time,  then 
$x  = 0 and  $ reduces  to  a function  of  time  only.  The  three  con- 
ditions X = Y = Z = 0 then  give, 

i 

(Y  - l)(c"  - q2  - 23>t)  q2  = 2q2q2 


2A2qqx  + (y  - DA^q2  = + 2qq+)  7 (4.18) 


A'  = 


v, 


The  functions  Aj_  are  given  directly  in  terms  of  q(X;t),  © ( X; t ) 
and  their  derivatives  by  equations  (4.8).  ^(t)  may  be  eliminated 

between  the  first  two  equations  of  (4.18)  to  yield  an  equation  for 
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fe 

i ■ 

\ 


q(X;t)  in  addition  to  the  third  equation,  A2  = A^A^,  This  equa- 
tion is 


2A2qq\  + - 1)Aif  qx  = A1  { 


qq i. 

t y - 1 dt 


2 2 

-q  v 


at  (-?r>  | C+-19) 


(4.18),  and  (4,19)  are  therefore  the  general  equations  for  the 

velocity  for  the  unsteady  simple  wave  when  the  "centre"  is  fixed. 

It  is  obtained  in  terms  of  any  parameter  X(x,y,t)  (for  which  VX 

does  not  vanish)  which  may  be  taken  to  be  the  velocity  q(x,y,t). 

The  original  functional  relation  q = q(X,t)  now  simply  reduces  to 

q = X,  so  that  q^  = 1,  q^  = qt  = q^  = q t - 0,  whilst  ©(X;t) 

becomes  Q(qr.t)  and  ©■»  = 6 , ©*  % = 6 , etc, 

^ X q5  XX  qq’ 

Equation  (4.19)  now  reduces  to 

q£lA1 


-Apq  + (y  • £ 


3 / “2 
(q  ) 


Y - 1 at  'i 


(4.20) 


where 


q2  = 1 + q2©q2 

^ 2 2 

1 + q e » _n 

A = 5 3_  -—-[©  + tan  q©  ] 

1 _27  3q  q 


q 

1 + a2©  2 


s 


(4.21) 


A2  ~ 


q2© 


3 r<^  • j-^„“1.  ^ t 

sr  ■r  I'**!!  uv  j 

3t  q 


- ©2  + :tt 


q 

etr 

q® 


A,  = t 


fl 

q© 


Equations  (4,18)  and  (4,20)  now  become  two  partial  differential 
equations  for  the  single  dependent  function  ©(q;t),  q now  being 
regarded  as  an  independent  variable.  Thus  it  appears  that  the 
flow  is  overdetermined.  The  above  simplification  of  the  general 


V* 


! 
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theory  may  be  obtained  directly  by  resorting  to  the  equations 
(2,1)  and  (2,2),  It  should  be  noted  however  that  this  direct 
approach  is  only  possible  when  the  A-lines  are  concurrent  through 
a fixed  point,  and  may  therefore  be  considered  as  radii  through 
the  origin  when  the  current  polar  coordinate  r (the  radial  dis- 
tance) is  of  no  consequence  to  the  velocity  £.  It  is  not  capable 
of  direct  generalization  to  reproduce  the  more  general  equations 
(4,10a)  for  which  there  exists  a time  varying  envelope.  In  view 

of  the  presence  of  three  equations  for  just  two  functions  ($,0) 

I 

or,  in  the  simpler  case,  .wo  equations  for  0,  thC  S 0 equations  will 
be  established  in  the  simpler  case,  by  a more  direct  method. 

Having  thus  verified  their  validity,  it  will  be  then  justifiable 
to  return  to  the  general  equations  (4,11),  (4,12),  and  (4,13). 


Direct  Procedure  of  the  Simplified  Case.  Let  u,v  be  the 
components  of  velocity  along  and  normal  to  the  radius  vector 
drawn  from  the  origin,  d the  radial  angle  and  = 0 - d the  angle 
q makes  with  the  radius  vector.  Then  u = q cos  ^i,  v = q sin  \Jj 
and  the  general  equation  for  the  velocity, 

at(ur  + £ + * vd)  = q(uqp  + J ) + 2q4  + 9 

reduces  to,  since  r is  of  no  consequence  to  q(d;t), 

r- ) 

ac'(u  + vd)  = vqqd  + 2rqq  + rq> 

that  is, 

(a2  - q2)qd  sin  ^ + a2q©0  cos  \Jj  = 2rq<5  + rep  (4.22) 

Now  since  q(d,t)  is  independent  of  r 
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reduces  to 


r ■ 

cp(r,cs$t)  = J q * dH  = q*I.  + $(t) 
cp  = O(t)  + rq  cos  \}> 


Also  the  condition  y*q  = i.e.,  v = u gives 

A—  o 


cot  \}j  = q© 


(4.23) 


from  which  by  differentiation  is  obtained 


q©  + q(©  + q9  ) = - cosec2  ij>© 

q qt  qq 

regarding  Q(d$t)  as  a function  of  (q;t)  so  that  © 
etc.  Therefore 


= © q + ©+, 


& 


o P 

q [q©  + (1  + cosec*^)©  ] + q©  + S cosec  ^ 

qq  q qt  t 


= 0 


Similarly 


q [q©  + (1  + coscc^)©  1 = cosecc'^ 

d qq  q 


The  expressions  for  the  time  derivatives  of  9 are 


m 

fe- : 


i 


m 


k 


'vtr. 

i 


9 = - rq©  sin  ^ 

t 


9 = + 


A(q;t)r  sin  *4> 


^ q9  + (1  + cosec2,|0© 

qq  q 

after  successive  differentiation  and  substitution  in  the  resulting 
expressions  of  q,  where 


A(q;t)  = (©.cosec2^’  + q©  , )2  - q(q©  + l+cosec2^*©  ) (©^^+©,  cotxb ) 
t qt  qq  q tt  t 

and  the  expression  for  the  acoustic  speed  becomes 
2 

= \ (c2  - q2)  - 9 + rq©  -sin  $ 


\<  'tv 
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When  these  expressions  are  inserted  in  (4.22)  and  the  coefficients 

p 

of  (1  : r : r ) equated  to  zero  there  results  the  three  relations 

fi 

(Y  - l)(c2  - q2  - 20  ) = 2q2sin2'|>  j 

^ ! 

(Y-l)©^  cosec2‘>5  + 2(qG^^+9^.  cosec2^)=q“^$^  (q©^+l+cosec2^»8^)  r 


A(a  ; t)  = 0 


(4.24) 


which  may  be  reconciled  with  relations  (4.10)  by  noting  that,  in 


(4.21),  since  d 


= 9 - = © - cot“^q©  = © •;  tan_1a©  h constant, 

q '4  ’ 


then  ^ 

2 , q©  ^ + ©^  cosec 

A = cosec  U>_  (j  = qt  t 

2 q cot  -vJj  2q 

9 q 


A = cosec~~^  d _ 
1 q cot  ^ q 


(1  + cosee^)©,- 


> (4.24) 


A = tan  ^(©  + © cot  ^) 

j;  "Cut 


A,  = © tan  ib 
4 t 


A(q;t)  = (A^  - A A )qV 
2 13  q 


Equations  (4.24)  therefore  reduce  to 

(Y  - l)(c2  - q2  - 2<P  Jcosec^  = 2q2 

Vs 

2 

(Y  - 1)A,,  cosec  ^ + 2qA0  = A^  $ 


(4.26) 


a2  - axa3  = 0 


These  coincide  with  (4„18)  since,  in  the  original  notation, 

2 2 2 2 
q,  = 0,  and  q - 1 + q © = cosec  ^ • 

X/  ua  n 


! 


' 1 


I 
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5.  The  General  Case.  The  following  analysis  deals  with 
the  case  when  the  parameter-  lines  do  not  necessarily  possess  a 
degenerate  envelope.  The  relations  (4.21)  are  transformed  to 
give,  in  terms  of  p = ©z;  q = ©t 

A1  ' £ {m  * P2)  - ezzK2z  a2  = 1 {q(l  + p2)  + eztj  e"z 

A3=HPq2+9«}  \ = -p- 

since  z = log  q;  q©  = ©z  = p.  Similarly  the  expressions  for  the 
functions  Bj,;  (4.8)  become, 


- eZB,  = i + 

1 P 

- B,  = 111*1 


-Z 

- e B = -44 
3 P 


- e 


”zBi  = 3. 


4 


\ 

ft 


where  q now  denotes  the  derivative  ©t,  and  is  not  to  be  confused 
with  the  total  fluid  speed  now  denoted  by  ez. 

The  angle  d which  the  parameter  lines  make  with  a fixed 
direction  © = 0 is  given,  apart  from  a constant  additive  angle, 
4y 


that  is 


d = 0 + tan-1©,.. 


d = © + tan_1p 


since  P , the  angle  at  vhich  the  streamlines  intersect  the  para- 
meter lines,  is  given  by 

cot  r = ©z  = p 


and  d = © - T, 
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Solution  of  the  Velocity  Equation.  Of  the  three  funda- 
mental simple  wave  equations  (4.11),  (4,12),  and  (4,13) * only  the 
last  two  involve  (£(q:t).  The  relation  (4,11)  involves  only  the 
velocity  and  is  independent  of  $>,  therefore  it  is  possible  to 
obtain  general  properties  of  the  flow  predicted  by  the  common 
equation  (4,11)  for  any  possible  form  of  *(q;t),  (4,11) 

p 

(A2  = A-^A^)  gives  for  the  equation  for  9(z;t), 

C©zt  + qd  + p2)]2  = [ett  + pq2]  * C ©2Z  + pd  + p2)l 


which  becomes 


pq2ezz+p<1+p‘2)s 


tt 


-2q  ( 1+p2  ) ©z  t-q2  ( 1+p2 ) =0  (5.D 


which  is  equation  of  the  Mange-Ampere  type.  It  is  non-linear 
but  the  second  order  terms  of  the  second  degree  occur  only  as  the 
Jacobian  9(p ,q)/9(z ,t ) = ~ ®zt  therefore  it  can  be  oealt 

with  by  considering  the  characteristics  which  are  given  by 
(cf,  Forsyth,  Differential  Equations:  p,  ’185)  • 

dp  + p(l  + p2)dz  + Xdt  = 0 

O 

dq  + \dz  + pq^dt  = 0 


where  X is  given  by  the  quadratic 

x2  - 2Xqd  + p2)  + q2U  + p2)2  = 0 

p 

that  is  X = q(l  + p ).  Thus  the  equation  for  9(z;t)  is  always 
parabolic,  the  roots  X being  real  and  coincident. 

The  coincident  characteristics  are  given  by 
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dp  + p(l  + p2)dz  + q(l  + p2)dt  = 0 


O 

dq  + q(l  + p£l)dz  + pq^dt  = 0 


wliich  mC'.y  be  ccinbx  ned  to  eliminate  dz  in  the  form 


qdp  - pdq  + q dt  = 0 


so  that,  if  u is  an  arbitrary  parameter  of  integration, 


t + £ = u 

q 


(5.2) 


Similarly,  elimination  of  dt  yields, 


pqdp  - (1  + p2)dq  - q(l  + p2)dz  = 0 


so  that  on  integration  in  terms  of  another  parameter  v, 

2 2z  2, 

q e = v ( 1 + p ) 


(5.3) 


It  follows  that  the  "Intermediate  Integral"  of  (5.1)  may  be  ob- 
tained by  taking  any  arbitrary  function  of  (u,v),  that  is 

F(u,v)  = 0 


which,  being  an  equation  (involving  p,q)  of  the  first  order  may 
be  further  integrated  for  9(;c,t),  if  desired  by  Charpit’s  general 
method  for  such  non-linear  first  order  partial  differential 
equations. 


The  general  intermediate  integral  may  be  written 


V = ? ?Z  9 

q'e'  = (1  + p2)f(t  + £) 


(5.4) 


where  f is  arbitrary.  However,  it  is  possible  to  evade  the 


i 


1 


Vi  n 
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general  method  by  solving  (5.2),  and  (5.3)  for  p;q  as  functions 
of  (z,t,U,v)  giving,  since  p = q(u  - t), 


q = - 


'/v 


J~2 z Z~ 172 

v 6 — ■ v\  — t ) 

v/v  • (u  — t ) 


- v(u  - t.)‘ 

where  actually  v = f(u),  an  arbitrary  function  of  u* 
Then  insertion  in 


gives 


d©  = pdz  + qdt 
de  _ (u  - t)dz  + dt 


/v  V>  - -a 


- v (u  - t) 


so  that 


0 - c = cos”--*-  -j  (u  - t)  v/v  e Z j 


therefore  ©(z,t)  is  given  by 


__  **  z 

cos  (Q-c)  = (u-t)'/ve 


or 


q cos  (0  - c)  = (u  - t)  /v 


where  q now  denotes  the  fluid  speed. 

If  c,  v are  regarded  as  arbitrary  functions  of  a single 
parameter  X(  = u)  then 


q cos  |©  - cp(  X)  j-  = \}j(  X)  ( X - t) 


(5.5) 


is  the  !l complete"  integral,  from  which  the  "general"  integral  is 
obtained  by  partial  differentiation  with  respect  to  X, 
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9’  (X)q  sin  |©  - cp(X)  j-  = (X  - t)^’  (X)  + ^jj(X)  (5.6) 

and  eliminating  X between  (5. 5)  > and  (5.6)  [definitely  chosen 
functions  cp(X),  ijj(X)  would  be  necessary  in  actually  carrying  out 
this  procedure]* 

It  follows  from  (5.5)  in  terms  of  z = log  q 

cos  [6  - cp (X )]  = i)j(X)(X  - t)e~Z  (5.7) 


therefore  partial  differentiation  with  respect  to  z gives,  after 
considering  the  corresponding  variation  of  X(z,0,t)  which  accord- 
ing to  (5.6)  vanishes,  giving 


therefore 


o 

z 


d 


reduces  to,  simply 


cot  [0  - cp  (X ) j 
0 - cot”1© 

& 

d = cp  (X) 


(5.8) 


(5.9) 


(choosing  principal  value  for  convenience),  d denotes  the 
inclination  of  any  local  parameter  line  (along  which  9,q  are 
instantaneously  constant)  to  6 = 0,  therefore  these  lines  will 
constitute  a system  of  parallel  lines  if  o(X)  reduces  to  a con- 
stant, Since  in  this  case  the  left  hand  side  of  (5.6)  vanishes 
and  the  right  hand  side  then  merely  implies  that  t is  an  arbitrary 
function  of  X,  then  X is  an  arbitrary  function  of  t so  that  (5.5) 
gives 

q cos(©  - ©0)  = F(t) 

F(t)  an  arbitrary  function  of  t and  ©0  a constant. 

The  combination  of  (5»5)»  and  (5.6)  is  not  the  most 
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general  solution  for  9(q;t);  a more  general  solution  may  be  ob- 
tained by  applying  Imschenetsky*  s method  by  inserting  in  (5.1) 

the  solution  +z 

e cos  (6  - c)  = b(a  - t) 

and  allowing  the  parameters  (a,b)  to  vary,  regarding  c as  a 
function  of  (a,b),  However  attention  v/ill  be  confined  to  the 
class  of  s orations  obtained  by  eliminating  X(q;t)  between 


£ 


q cos  (9  - X)  = «p(X)T4»(X)  - t]  (5.10) 

q sin (9  - X)  = cp'  (X)  Ji|>(X)  - 1 1 + <p(X)\|>»  (X)  (5.11) 

which  are  obtained  from  (5*7) j and  (5.8)  after  a slight  change 
of  notation.  With  this  notation,  X becomes  equal  tod  ( equal  the 
inclination  of  any  parameter  line  to  9 = 0) , and  is  defined 
implicitly  in  terms  of  the  arbitrary  functions  cp  (X ) - ^ (X ) by 
means  of 

q2  = (cp2  + cp^)(\|j  - t)2  + 2cpcp^i|)^(\{j  - t)  + 92^2  (5.12) 

obtained  from  (5.10),  and  (5.11)  by  elimination  of  9.  It  is 
always  possible  to  rewrite  the  system  (5.5)>  (5*6)  in  the  form 
(5*10),  and  (5.11)  provided  that,  in  the  original  system  the 
function  cp(X)  [occuring  in  (5.5)  > and  (5.6)]  does  not  reduce  to 
a constant. 

Since  0 

q9q  = cot  a = Qz 

q©£  = cp (X.)  cosec  a 

a = 9 - X 

i i 

X(q;t)  = 9 - cot  q9  I 

q i 


(5.13) 


i 
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q9  = - (q  cot  a - X )cosec^a 
^ zz  H q 


57 


q20 


tt 


9^cosec2a  cot  a + q cosec  0(9^  + 9 cot  a)X^_ 

^2, 


(5.1*0 


q©  = - cosec 9 cosec  a - qX 
Zt  ^ 


1 

tj 


The  partial  derivatives  of  X(q$t)  may  be  obtained  from  (5*12), 


giving 


LAX  = K 

v 


LAXq  = q 


where 


L = - t)  + qn})^ 


K = (cp2  + cp^)(i|)  - t)  + 99^ 


A = (9  + 9^)(^  - t)  + 2w.^x  + qnj)^ 
so  that  (5.1*0  becomes,  since, 

q cos  a = 9(tJi  - t) 
q sin  a = 9^(^  - t)  + = L 


f 


A 


then 


X 

o 


n zz 
'*1 

IT  + 

9 - 

tr-  0 

n 

j^ZZ 

= q2[l 

- <POI» 

- t)A]  ! 

1 

| 

< ^zt 

i 

] 

= q2[K 

- 9A] 

L 

S' 

t 

! L3i9tt 

- K2  - 

cp^  (4>  - 

t)A 

1) 

(5.15) 


(5.16) 


(5.17) 


(5.18) 
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& 


£ ■'>* 


The  expressions  (*+.21)  for  A^, 


2 


qA^_  = cosec  a tan  a X^ 

p 

qAg  = cosec  a tan  a X 


qA.  = tan  a cosec  a ( cp^  + cp  cot  cOX^  j 

I 

{ 

qA^  - cp  tan  a cosec  a J 

become,  with  the  aid  of  (5.1V).  and  (5.15) , 

o o 0 

irA<p(4i  - t)A  = q^  i 


*-• 


therefore 


L2Acp(4>  - t)A2  = qK  l;> 

L2Acp(\}>  - t)A^  = K2  1 
(t|>  - t)\  = 1 J 


A1  - fS.  = tl  = = 


c?  qK  K2  L2Acp  L2Acp  (i|)  - t) 
It  is  clear  from  (5.19)  and  (5.21)  that, 


^3.  _ k2  _ At 

at  - t - r? 


§1 

and  that,  from 

(5.15), 

F 

& 

xA 

xv. 

Q 

X 

II 

O 

W.  • 

q 

* t 

W& 

ft 

r 

% 

ttCZ? 

sr 

therefore 

Also; 

A-.X2 

3 q 

- A-,X2  = 0 
1 t 

i 


T,nOq  = 9(.|.  - t) 

LQt  =9 

q2  = 1 + q29q2  = -1. 


(5.19) 


(5.20) 


(5.21) 


(5.22) 


(5.23) 


i 
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Av  • 


Relations  Between  If  z = log  q then 

2 


q2A  = 1 + Q2  + 

^ 1 z © 

i z 

) a3  = »t  + 5^ 
z 


qAj  = V1  i SU-g^  1 


©, 


©f 

A = 

4 . 


n e 

- qBi  = 1 + _zz 

J z 

] - 1b.  = !it 

i,  q 3 ®z 


. B = °t  +-  “zt 


~i  B = Si 

q 4 ©. 


J 


Hence 


j\  = - V 


b3  “ «<et  - V 


B2  = «<Vt  - V 

\ = " C1A4 


(5*24) 


I- 

* 


fe 

1 

9^ 


at 


There  Exists  the  Following  Identities,  If  we  eliminate 
(0  ,0 x)  between  equations  (5.24),  using 


then 
that  is 


(q©2)(qQ2)  = (q©q©t)2 
(qA-j_  + B1)(qA^  + B,  ) = (qAg  + B2)‘ 


(B1B3  - b|)  + q2(AxA3  - A2)  + q(AxB3  + A^  - 2A2B2)  = 0 

Also,  directly, 


fj 

9 


iv 

w 


: i 


i(A;LB3+A3Br2A2B2)=  -2(A1A3-A2)  + (A10t-2A20q©t+A3e^) 


i-VBlB3-B2)  = <A1A3  - Ai>  - ^AlQt-2A2©q©t+A3eq) 


(5*25) 


q 

j 
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The  Linear  Equation  for  <Kq;t).  The  two  Monge -Ampere 
equations  (*+.12),  and  (4.13)  are  both  equations  for  $(q;t)  in 
which  the  coefficients  A^ , are  functions  of  (q,t)  and  Q(q;t), 
the  solution  of  (5.1)  as  expressed  by  the  system  (5.10),  and 
(5.H).  The  final  expressions  for  the  coefficients  Aj_  are  given 
in  equation  (5.21)  in  terms  of  arbitrary  functions  <p(>,)>  i|j(\) 
the  variable  \(q;t)  (given  by  (5.12))  and  the  functions  L(X;t), 
A(X;t),  K(X;t)  (given  by  (5.16)).  The  equation  (4,12)  is  a 
general  second  degree  Monge-Ampere  equation,  but  (4,13)  is  of 
the  first  degree  in  ^(qjt).  It  is  preferable  therefore  to  study, 
first,  solutions  of  (4.13)  end  then  return  to  (4,12)  to  obtain 
any  possible  solutions  common  to  both  these  equations. 

The  equation  (4.13)  may  be  written 

2 

A <I>  -2A0$  +A,  ■ + (A,B.  +A-B--2A,B0)<I>  =2qA0+(y-l)A,  q 

3 qq  2 qt  1 tt  ' 3 1 1 3 2 2 q 4 2 4^-q 

p 

or,  in  virtue  of  (5.25),  since  A-^A^  = A^, 

A3hq-2A2*qt+Al®tt+Cl0q<A3eq-2A2eq9t+Aie?)=2o-A2+(r-l)A^  (5.26a) 

and,  using  (5.21),  and  (5.2j) 

K2<5qq-2qK$qt+q2(Dtt+q(Pq(Keq-qQt)2  = q2  j2K+(y-l)cpA  j 
But  from  (5.23) 

KG  - q0t  = JL  [KOI)  - t)  - q2]  = - 1-lX. 

h u qL  q 

after  simplification  with  (5.16). 


- r ' T^-ta4 
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The  final  equation  for  0 is  therefore 

K2$  - 2qK$  , + q23>  + 9^  1 9 = q2  { 2K  + (y  - 1)9^1  (5.26b) 

Qt  bt/  ^ Q 4.  ^ ' 

This  equation  is  parabolic , the  coincident  characteristics  being 


given  by 


that  is 


q2dq2  + 2qKdqdt  + K2dt2  = 0 


qdq  + Kdt  = 0 


(5.26c) 


which  in  virtue  of  (5.1 5)  may  be  written 


X dq  + X. dt  = 0 
q t 

Therefore  the  characteristics  coincide  v/ith  the  curves 


X(q-,  t)  = constant 


(6.27) 


Form  (5.12),  these  X-curves  constitute  the  one  parameter  family 
of  hyperbolas 

(cp2  + <p2)(il>  - t)2  + 2cocp^x(il)  - t)  + 924'J  = q2 

in  the  (q;t)  plane.,  The  subsidiary  system  of  the  equation 
(5.26)  is 

K2dPdt  + q2dQdq  = [q2(2K  + y~-  1 cpA)  - cfH^Pq"1  ]dqdt 
where  P r-  $ ; Q = © • V/ith  the  use  of  (5.26c),  this  system  reduces 

'l  V 

t0  K 9^^  

- ({*)  dP  + dQ  + — Pdt  - (2K  + y-  lcph  )dt  (5.28) 

q qJ 

Since  X(q?t)  = const,  are  the  characteristic  curves,  it  is  possible 
to  integrate  the  system  (5.28)  regarding  X,  9(X),  \j)(X)  as  con- 


r- -v3*  r- 


stants 
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I 

i" 


V 

K> 

& 


K' 

& 


: it 
*5  . 


cd 


i l-v 


For  this  purpose  (5.16),  and  (5.12)  give 

2 


_ k 

dt  q 


Then 


»a  A <5>  ■ I - K & 


therefore 


,3  d /K' 


4.  2 


q'  at  (5>  = K - q‘'<(f  +<PX>  = - <P  \ 

Insertion  into  (5.28)  enables  the  integration  to  be  performed 


giving 


KP 


[2K(X,t)  + y - 1 cp(X)A  (X  ,t)]dt  + 0 


where  the  indefinite  integration  ic  to  be  performed  v/ith  X con- 
stant, and  C is  a constant  of  integration.  The  intermediate 
integral  of  (5.26)  is  obtained  by  making  C an  arbitrary  function 


of  X.  therefore  v 

4>t  - q % = C(X)  + F(X,t) 


(5.29) 


where 


-F(X,t)  = J Ky  + l)?1"  + 2cpJ  -r  (y  - l)<pcpxx  1 - t)2 


+ cp(2Y<p^i]^  + Y “ 1 'I1-  t)  (5.30) 


This  intermediate  integral  is  linear  in  O and  has  the  subsidiary 
system 


d<£> 


dt  _ 9Ciq  

1 “ K C(X)  + F(Xst) 


in  which  X(q$t)  is  still  given  by  (5.12),  The  first  of  the 
equations  implies 


qaq  + Kdt  = C 


which  (cf,  equation  (5.26c))  defines  the  original  characteristics 


I 

I 
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already  encountered,  and  may  be  taken  as  a first  integral  of 

(5.29) 

X(q;t)  = constant 

and  therefore  this  condition  may  be  employed  in  conjunction  with 

= C(X)  + F(X;t) 

to  yield  the  second  integral  of  (5.29)  (by  partial  Integration 
with  respect  to  time  with  X constant)  in  terms  of  a constant  D. 

r 

$(q$t)  = - C(\)(4)  - t)  + | F(X',t)dt  + D 

u 

The  general  integral  of  (5.29),  and  therefore  the  complete  inte- 
gral of  (5.26),  may  be  obtained  by  making  D an  arbitrary  function 
of  X,  so  that 

<&(q;t)  = D(X)  - C(X)0|>  - t)  + E(X)(t|>  - t)2  + G(X)0l>  - t)3  (5.3D 
where  C(X),  D(X)  are  arbitrary,  and 

G(X)  = 1 [(y  + l)cp2  + 2cpg  + (Y  - l)99xx] 

E(X)  = i 9[2y9^)x  + (Y  - 

The  i-tonge-Ampere  hquation  for  $(q;t).  The  relations 
(5.10),  (5.11),  and  (5»31)  define  0(q;t);  $(q$t)  in  terms  of  the 
arbitrary  functions  cp(X),  ^(X),  C(X),  D(X)  of  the  parameter 
X( q 5 1 ) given  implicitly  in  terms  of  9(X),  \j)(X)  by  means  of 

q2  = (cp2  + 92)(i'j  - t)2  + 2cpc|y})x(4>  - t)  + 9^\  (5.32) 

It  is  necessary  to  find,  therefore,  if  it  is  possible,  by  suitable 
choice  of  these  four  arbitrary  functions,  to  satisfy  the  additional 
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c-H- 


relation.  (4.12)  between  <3>,  0.  This  relation 

(V  + Bl4’,K<Stt  + B3V  " %t  + B2<t>q  + 2q>(*qt  + B2V 


- q + (Y 


may  be  written 


- - \ - Vq}(1  + ^ = 


0 


- 2q<»qt  + ®q  (B3*qt  - 2B2  • *qt  + Vtt5  - 2qB-* 


qq  tt  qt 


2 q 


(B^  - B|)^  - q2  + (Y  - 1) 


2 

- - Vq)  (1  + ^ = 0 


Like  equation  (5«1)  this  is  a Monge-Ampere  equation  for  <P C q 5 1 ) 
in  which  the  coefficients  are  functions  of  q and  9(q;t)  as  well 
as  since  the  second  order  terms  of  the  second  degree 

occur  only  through  the  presence  of  the  Jacobian  3 (^q;<J>t)/d(q;t) 

r\ 

= “ qt*  Th:*-S  Monge-Ampere  equation  may  be  reduced 

slightly  by  using  the  relation  (5*1)  (A^A^  = for  and 

the  relation  (5.26)  for  <&(q;t). 

Substitution  of  the  expressions  (5.24-)  for  enables 
the  equation  for  $ tc  be  written 

- $2,  - 2q<5  . + q$  (©20  - 29  9,<&  + Q2<5  , ) 

44  tt  qt  qt  q t qq  q t qt  q tt 


- «YA3%  ' 2A2*qt  + AAt>  - 2q2(eq®t  " V*c 

[2 


2-2 


- " 2A2Vt  + A39q>  - + <Y  - 1)  {— s-3- 

- «t  + SV’q]  (1  + = 0 


Now  using  directly  the  linear  equation  for  $ (equation  (5- 26a)) 
this  may  be  reduced  to 
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fT\  A} 

“qq  'tt  "qt 


$2  - 2qO  + q<J>  (92$  - 20  6 4 + ©2  $ ) 

qt  qt  q t qq  q t qt  q tt 


- q2(l  + 26q9tg>q)  + (Y  - 1)  [S— - ®tJ(l  + q2©q)  = 0 

(5.33) 


that  is 


(4qq  + t‘eq4q>  (*tt  + 4®tV  ' (V  + qVt°q  + ^ 

* (r  - i)  {sif-ai  - \)  a + ^ = 0 

(5.3*0 

which  does  not  contain,  explicitly,  any  of  the  functions  (A^,B^) 

nor  the  second  derivatives  (9  . 0 . G ) but  only  (0  , 9,  )• 

qq7  qt7  tt  q t 

The  equation  (5«33)  (or  its  equivalent  (5»3'-0)  is  not 
parabolic.  The  characteristics  will  be  obtained  by  first  finding 
the  roots  of  the  equation  in  p,, 

l 2 2 ) 

p2  - 2q(l  + 9q9t$q)p  + q2 (1  + 9q9t<5q)2  = (Y-l)|— -®tJ(l+q292 ) 

(again  cf,  Forsyth,  Differential  Equations,  p.  485)  that  is,  by 
the  roots  p., 

* 2 2 

[u  - q(l  + 9 9.  <6  ) I2  = (v  - 1)  -!  - — "—2 <P  l (1  + a2Q2) 

L-  H q t q - • 12  t j q- 

Now  the  equation  (5»34)  will  be  parabolic  only  if 

- 1 


or 


<£  = J (c  - q ) 

®(q;t)  = | (c2  - q2)t  + F(q) 


|& 


but  resubstitution  into  (5.3*0  gives,  since  $ = $ , + q = 0, 

XX  QX  * 

0 0 0=0 
qq  q t 
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9,  “ 0 is  a degenerate  (steady)  flow  whilst  O . <X>  cannot  vanish 
t q ' qq 

and  be  simultaneously  compatible  with  = l/2(c2  - a2). 

In  the  general  case  it  would  be  necessary  to  insert  In 
(5.3*+)  the  previous  solutions  for  ©(q;t),  3>(q,t)  in  terms  of  the; 
functions  <p(\),  ^(M,  C(\)  and  D(\)  which  must  therefore  be 
chosen  to  satisfy  together  the  resulting  equation*  This  procedure 
would  be  very  complicated,  however-  owing  to  the  expressions 
for  ©n , ©^  and  the  derivatives  of  \(q;t)  invoiced  during  the 
substitution.  Certain  important  consequences  may  be  obtained 
from  a direct  consideration  of  the  general  solutions  for  © and 
<3>  without  any  further  regard  to  the  conditions  to  be  satisfied 
by  the  four  functions  cp,  i}),  C,  D,  Therefore  we  shall  not  carry 
the  general  theory  any  further  but  merely  state  that  any  solution 
of  a generalized  simple  wave  must  necessarily  be  of  the  form 
predicted  by  equations  (5.10),  and  (5«3D. 


Concurrent  Parameter  Lines.  If  the  lines  <J>n  + r • q = 0 

" u — q 

are  concurrent  through  a fixed  point,  which  without  any  restric- 
tion may  be  taken  to  coincide  with  the  origin  r = 0,  then  0 
reduces  to  a function  of  time  only  and  equations  (5.26b)  and 
(5.3*0  reduce  to 


%A,  = 2K  + (y  - 1)?£ 

r_2  2 


rtt 

2 


■ (r- 


where,  from  (5.23) , a2  = qVL2,  therefore 

q 

p -2  ,2 


*t“ 


Y - 1 


! 
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Inserting  in  these  expressions  for  0^.;  the  values  of  L,  K, 

£ in  terms  of  <p(X),  \j)(X)  as  given  by  (5»16)  gives 

+ 1)9'  + 2<PX  + <Y  ~ DWxxlC*  “ t)+2Ycp(px4>x+(Y-i)'p2\x 
2«p>  r i 

®t*  = riK(^  - t}  + j 

Now  assuming  that  X does  not  reduce  to  a function  of  t only 
(otherwise  (5.12)  will  give  the  contradiction  that  q is  a function 
of  t)  the  above  simultaneous  expressions  for  which  is  a 

function  of  t only,  must  simply  be  of  the  form  (a  + bt),  a,b 
being  constants  independent  of  X or  t. 

These  give  the  four  conditions,  in  terms  of  constants 

(a  jP ) , 2 2 

(y  + D?  + 2<px  + (y  - 1)9<PXX  = o 


. + Cy  - D<j4xx  - P 

2cpx  = (Y  - 1)0 

299x^x  = (Y  “ 1)P 

These  equations  are  clearly  incompatible.  The  third  gives  at  once 

?xx  = 0 and  

9(X)  = ± X + cpQ 

which  enables  the  first  tp  be  reduced  to 

n f M - V ^|2 

(Y  + 1)  9n  ± X -V  + (y  - 1)«  = a 

b u V ^ J 

which  is  an  impossible  relation,  even  if  a = 0 then  c?Q  = cp  = 0, 

It  is  concluded  that  there  is  no  unsteady  generalization  of  the 
simple  wave  for  which  conditions  are  instantaneously  constant 
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along  a series  of  straight  lines  radiating  from  a fixed  point. 
However,  this  result  says  nothing  about  the  possibility  of  lines 
concurrent  through  a moving  point. 

Oscillating  Flews.  The  principal  possible  application 
of  the  preceding  solution  would  be  found  if  it  were  possible  to 
choose  the  four  functions  cp,  0,  D in  such  a way  that  the  time 

variation  of  the  deflection  function  ©(q;t)  was  harmonic.  Thus 
consider  a flat  plate  whose  equation  is  written  in  terms  of  the 
parameters  a(t),  n(t)  (the  unit  normal  vector)  in  the  form 

a + £ • n = 0 (5.35) 

then  the  boundary  condition  to  be  satisfied  on  this  plate  is 


[at + a 


V)' 


r * n)  = 0 


giving 


+ q • n *r  r • n 


n = 0 


(5.36) 


In  addition  there  is  the  field  equation,  valid  everywhere, 


$L  + r • q =0 

4 " “1 


(5.37) 


In  the  first  place  it  is  noted  that  the  line  (5.35)  cannot  ) 

< 

instantaneously  coincide  with  a parameter  line  of  the  family 

(5*37),  except  for  the  case  when  (5.35)  simply  represents  a 

translational  oscillation  of  the  line  with  no  pitching  (i.e,  n 

- a constant  vector).  For  otherwise  (q,©)  would  not  be  constant 

along  such  a line,  as  required  by  the  boundary  condition  (5»36) 

* 

when  n = 0.  Thus  in  the  general  case  it  is  possible  to  eliminate 
r from  the  equations  (5.35)?  (5.36),  and  (5. 37)  to  give 


. a- 


All-111 


69 


= ~ a^2.  ) (a  f n ‘ q)(n  2.  ) (5A0) 

q -A-  A q a “A^q 


If  n = (cos  a,  sin  a);  a - a(t)  = ut  and 


q = (q  cos  9,  q sin  9) 

then 

q = (cos  9 - q9  sin  9,  sin  © + q9  cos  9) 

**q  q ’ q 


Therefore 

|n  q | = sin  (9  - ut)  + q9  cos  (9  - ut) 
A“q  q 


n • q = q cos (9  - ut) 


n^£  = - uf cos (. P - ort)  - q9q  sin(Q  - ut)] 


the  required  condition  (5«i+0)  becomes 


u£>  = ua[cos(9  - ut)  - q©  sin(9  - ut)] 

- [a  + q cos(0  - ut)]-[sin(9  - ut)  + q9^cos(9  - ut)] 


Now 


_ <p(\}>  - t ) _ 


= cot(9  - X) 


q cos(9  - X)  = cp(ij;  - t) 


q sin(9  - X)  = L 


and  if  a = aQ  sin  ut  then 


t J$>  ~ - ua  [sin(9  - 2ut)  + q6  cos (9  - 2ut)] 


- q cos (9  - ut)[sin(9  - ut)  + q9^cos(9  - ut ) ] 
finally,  on  inserting  q©  = cot (9  - M , 

'■i 


'.Sk^'W 
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sin(9  - 2ut)  + ad  cos(Q  2ut)  — — ■=— - — — — — - ^ cos(X  - 2ut) 
‘ sin(O-X)  L 

sin(6  - ut)  4-  a©  cos(Q  - cot)  = cos(X  - ut) 
q 1 

q cos(9-  ut)  = 9(^i  - t)cos(\  - ut)  - L sin(X  - ut) 
therefore  since  ^ = O^(^) 


- urtyr1 
giving 


ua 


cos(X  - 2ut)  + cos(X  - ut)  J 9(i(j  - t)ccs(X  - ut) 

' U 

- L sin(X  - ut)J- 


- u$^A  ^ = [ uaQcos  X - ^ L sin  2X  + ^ cpC^  - t)cos2\]  cos  2ut 

+ [uaQsin  X + ^ L cos  2X  + i-  9(^1  - t)sin  2X]sin  2wt 

+ I 9(9  - t)  (5.41) 

Now  having  regard  that  the  expressions  (5*31)  for  <p  and  (5.16) 
for  A are  simply  (finite)  polynomials  in  the  time  with  coeffi- 
cients functions  of  X , a representation  of  the  type  (5.41)  for 
$ can  only  be  possible  if  both  the  factors  of  the  sin  2ut, 
cos  2ut  terms  vanish  identically  in  X.  It  is  clear  that  this 
is  never  possible  except  in  the  trivial  case  aQ  = 9=0,  giving 
also  q = L = A = 0. 

In  the  exceptional  case  when  n = 0 then  (5-35) 5 (5*37) 
give,  since  they  must  represent  instantaneously  coincident  straight 
lines , 

aq  = 

— u 

■ q 

where  n is  a fixed  vector.  Therefore  q is  parallel  to  a first 

-q 

vector  and  the  system  of  lines  of  constant  velocity,  lines  (5*37)  > 


'•  jr*  ***  w •rarypyr'-rr*  ~ 


All-111 


71 


constitute  a system  or  lines  perpendicular  to  the  fixed  direction 
of  ft.  The  required  expression  for  $(q;t)  then  reduces  to,  if 
f(t)  is  an  arbitrary  function, 

•*»*>  ■ - 4 • + f(t) 

Possilbe  flows  falling  into  this  pattern  are  dealt  with  in  the 
next  section. 


t 

t 


6,  Parallel  Parameter  Lines.  This  section  considers 
the  class  of  flows  which  fall  into  the  analytic  form  £(\jt)  but 
which  vcre  automatically  excluded  in  the  analysis  following 
equation  (5,10),  Such  flows  will  be  shown  to  be  characterized 
by  the  velocity  being  constant  along  a series  of  parallel  straight 
lines  parallel  to  a fixed  direction,  and  the  component  of  veloc- 
ity in  this  direction  is  instantaneously  constant  throughout  the 
flow  field  and  is  a function  of  time  only.  These  represent  a very 
restricted  class  of  flows  and  it  is  concluded  that  the  previous 
analysis  is  always  applicable,  except  in  the  present  cases  for 
which  the  parameter  lines  degenerate  into  a system  of  parallel 
lines.  For  completeness  the  present  analysis  shows  the  procedure 
to  be  adopted  in  this  special  case, but  details  of  the  flows  are 
not  discussed. 

The  analysis  following  (5*10)  was  performed  on  the 
assumption  that  the  arbitrary  function  <p(\)  occurring  in  (5,7) 
did  not  reduce  to  a constant.  Otherwise  it  is  not  possible  to 
take  <p(M  as  the  current  parameter  to  replace  X and  rewrite 
(5*7)  in  the  form  of  (5,10),  If  cp(X)  is  a constant,  which  may 
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clearly  be  taken  as  aero,  then  (5*7)  becomes 

q cos  9 = (X  - t)\J>(\)  (6.1) 

from  which  a more  general  solution  for  9(q;t)  is  obtained  by 
partial  dif feriatiation  with  respect  to  X, 

0 = i|)(X)  + (X  - t)4»*  (X)  (6.2) 

and  eliminating  X between  (6.1)  and  (6.2).  Equation  (6.2)  implies 
simply  that,  since  \Ji(X)  is  arbitrary,  X is  an  arbitrary  function 
of  time.  Therefore  (6.1)  gives 


q cos  9 = qQ(t ) 

qQ(t)  being  an  arbitrary  function  of  t.  It  follows  by 
differentiation  that, 

Q„  = q9  = cot  9,  q9,  = - q (t)cosec  9 

p 

9 „ = - cosec  9 cot  9 
zz 

q9  . = q* (t )cosec^Q 
zt  o 

2 tt  | 2 2 

u 8^  = - qqo(t)cosec  9 - qQ(t)  cosec  9 cot  9 
The  expressions  for  A..  , become 


A,  = A-  =0,  A-,  = 


- qQ(t) 


qo(t)  ' 


A4  = 


qc(t) 


Since  = A2  = 0 it  follows  from  (4.8)  that 


9 + t,an~'Lq9  = 

'■i 


(6.3) 

direct 


(6.4) 


constant 
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and  the  constant  is  not  even  a function  of  the  time.  Therefore 


the  isovels  constitute,  at  all  times,  a system  of  parallel  lines 
parallel  to  a fixed  direction.  These  flows  could  then  clearly 
be  obtained  by  a more  direct  method  in  which  the  velocity  (but 
not  the  potential  or  the  acoustic  speed)  is  independent  of  one 
coordinate  x(or  y). 

When  the  above  values  for  are  inserted  in  the  linear 
equation  for  <&;  equation  (5«26a)  it  reduces  to 


®qq  + ~ , 2 2 \ ^q 

q(q  - q0)  4 


_ vY  - l)q 
2 T 

q - q~ 


that  is 


k v«a  - 4 


9 n. — a 

vq  - qD 


therefore 


q0(t> 

ii , , . 


q„(t) 

~n 

qD(t) 


0q  - (T  - 1) 


qqD(t) 

— n 

qo(t) 


qG(t) 

rs  2 

J a - a 


finally 


®(qjt)  = (y 


q2q‘(t)  r?  2 

- 1)  — ^ + G(t)  \lqz  - qi  + H(t ) 


2q0(t) 


(6.5) 


in  terms  of  functions  of  integration  ©(t);  H(t). 

These  functions  G(t)°,  K(t)  are  not  arbitrary  however, 
since  (6.5)  must  satisfy,  if  possible,  the  non-linear  equation 
(5.3l+)»  To  express  $(q|t)  in  terms  of  (yjt)  resort  must  be  made 
to  the  fundamental  field  equation,  which  relates  the  two  systems, 


Since 


+ r • q =0 

“q 


= (q  cos  0.  q sin  a)  = [q  (t);  q sin  ©] 
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then 

q = 

“q 

that  is 

n = 

Jq 

in  virtue  of  (6J1-), 
simply, 


[C  ; sin  9 + q0  cos  9] 

q 

[0  ; cosec  9] 

Therefore  the  field  equation  reduces  to, 
y = - $>  sin  9 

q 


or,  substituting  for  9(q;t)  from  (6.3), 

= „ 1 


v 


! 2 2 
^ - *0  * 0 


(6.6) 


q w o q 

Equations  (6.6).,  and  (6,5)  may  be  used  to  obtain  the  function 
q(y;t)  explicitly;  since 


12  2 


therefore 


y = - (Y  - 1)  \;q  - q*  • Hr—  - 0(t) 

>2 


qo(t) 


2 2 
i = cT  + 


u 

qjtr 


(y  - D2qo(t)'“ 


[y  + G(t)] 


(6.7) 


If  this  expression  for  q = q(y;t)  is  inserted  in  (6.5) 
an  expression  of  O as  a function  of  (y;t)  will  be  obtained.  This 
representation  is  of  the  form 


®(y;t)  - - h y2A(t)  + B(t) 


(6,8) 


where  A(t);  B(t)  [being  given  in  terms  of  qQ(t);  G(t);  H(t)] 
are  not  arbitrary  but  must  satisfy  the  equivalent  of  equation 
(5*?>4),  when  the  q-derivatives  of  0 are  transformed  into  ones 
with  respect  to  y. 

Returning  to  the  original  equation  for  q(y;t); 


t sr* 


All-Ili 


p o Q O • • 

(a  - u^u-x  - 2uvuy  + - v )vy  = 2(uu  + vv)  +9 


this  reduces  to,  since  u = u(t)  [=  qQ(t)j  and  v = v(y;t), 


2 2 • • „ 
(a‘  - v )vv  = 2(uu  + w)  +9 

v 


(6.9) 


Tt  should  be  recalled  that  the  condition  for  irrotational  flow 


v = u , is  automatically  satisfied  by  this  system  Fu(t),  v(y;t)] 
x y 


Also 


a"  _ 1 (2  2,  . 

y-1  2 C ~9)“9 


therefore, 


9 - <t>  + xu(t)  + yv(y;t) 


<5  =0;  + vv  •■=  0 

x y y 


9 = <E>t  + xu  + yvt 


* = *tt  + XUtt  + ^tt 


Substitution  into  (6.9)  gives 


(Y  - 1)  C ^(°2  - u2  - v2)  - - xut  - yvt]v 


- V vy  = 2(uut  + wt)  + 0tt  + xutt  + yvtt(6.10) 


u,  v,  <*>  are  independent  of  x,  therefore 


(Y  - l)u+vv  + u •=  0 


Put  u contains  t only,  therefore, 


so  that 


V = - s v 

7 (v  - 1 )u+ 


- Aft) , say 
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v(y;t ) - yA(t ) + D(t) 


u(t) 


exp  [ - 


(r 


i) 


^A(t  )dt  ]dt 


(6.11) 

(6.12) 


When  (6.8),  and  (6.11)  are  inserted  in  (6.10),  there  results  a 
quadratic  expression  in  y and  the  coefficients,  which  are  func- 
tions of  t only,  must  be  chosen  so  that  the  equation  vanishes 
identically.  This  will  yield  three  equations  for  A(t);  B(t); 
D(t).  [The  functions  A(t)  occurring  in  (6.8),  and  (6.11)  are 
identified  in  virtue  of  5 ■<-  yy  = 0,  ] These  equations  are 

i y 

Att  + (Y  + 3)AAt  + (y  + DA3  = 0 (6.13) 

Dtt  + (Y  + l)A(Dt  + AD)  + 2DAt  = 0 (6.14) 

(Y  - l)A[i  (c2  - u2  - D2)  + bJ  = AD2  + 2DDj_  - B,  . + 2uu.  (6.1?) 
t±  ° u Z Z Z 


The  first  equation  (6.13)  determines  A(t)  without  refer- 
ence to  the  remaining  functions  D(t),  B(t).  Let  X - 2/(y  + 1) 
and  rewrite  (6,13), 

\A  + 2 (X  + 1) AA  + 2A3  = 0 (6.16) 

which  may  he  factorized  into  either  of  the  forms 

(2A  + \ A)(A2  + |A)  - o (6.17) 

(2A  + ^)(a2  + x ~)  = ° (6.18) 

To  solve  write  p = dA/dtj  p dp/d.A  - d A/dt  and  w = AVp  so  that 


2AdA  = pdw  + wdp 
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then  (6,16)  becomes 


(to  + i)(u  + X)dp  + (!  + «+  X)du  = 0 


giving,  on  integration  in  terms  of  a constant  jx,  and  resubstitu- 
2 

ting  a = pw;  p = At; 


(A2  + At)^  = p,  (A2  + \At) 


(6,19) 


Since  k = 5/6  for  y = 7/5?  attention  would  have  to  be  confined 
to  the  special  values  of  u = 0,  ao  , since  it  is  not  clear  how  to 
integrate  again  the  first  order  equation  (6,19),  However  a 
parametric  solution  of  (6,16),  in  which  A and  t are  expressed  in 
terms  of  a common  parameter  z,  may  be  obtained  which  avoids 
resorting  to  (6.19), 

Let  A = then  At  A2  = Ytt/y  50  that  (6.17)  yields 


<2*t  - x'!'  = 


♦ 


that  is 


+4 


X -_£tt  + (2  - \)  It  _ 


o 


s. 

J 


tt  T 

therefore,  if  jx,  \Jio  are  constants, 

r iShrll  2JX-JJ.  ) 

• 2 2 ) i X i X ■ 

= » 1+  - +0 

When  y = 7/5  J X = 5/6  therefore 

,,2  _ 2, .-2/5  .-2/5. 

- I-1  - '!>0  ) 

5 

If  ^ = sin  9 then 
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6/5  P 5 

lit  = 5^0  j sin  ©dS 

giving 

6/^  U.  2 

3)it  = - i|)  'cos  0(3  ces^C  - 10  cos  © + 15) 

(6.20) 

or 

3|it  = - i^^zHz1*  - 10z2  + 15) 

(6.21) 

where  z 

= cos  0,  = t|>  (1  - z2)5/2. 

It  remains  to  express  A also  in  terms  of  ©, 

A = = 5 cot  ©0. 

4,  t 

where,  from  (6.20), 

.,,6/5  . 5^ 

Pv0  sm  ©©t  = v 

so  that  „6/4  5 

A(0)  = ' 'cos  0 cosec  © 


Therefore  the  sola  tlOii  Ox  \U«xuy  iiiciy  iic  &x  i/CUigCu.  3L* 


- — — *1  ^ J J »«*  4-  a h ^ . 


A = a cos  0 cosec  0 


(6.22) 


at  = - 5 cos  0(1  + 4 cos^G  - 4 cos20) 

5 3 <j 

where  a(=p/t)f  6^)  is  an  arbitrary  constant. 

When  these  results  are  inserted  in  (6.14)  in  which  D is 
now  regarded,  as  a function  of  G so  that 


5 

5 sin  © = aDg 

2 5 sin"©  = ac'(sin  - 5 cos  0 D0) 

then 

n . + 7 cot  0DO  - 10  D = 0 
Let  z = coo  0 then  the  equation  for  D becomes 

(1  - z2)D  - 3z  D - 10  D = C 

ZZ  z 


(6.23) 


(6.24) 
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Finally  (6.24-)  reduces  to  the  hypergeometric  equation 

p 

by  the  substitution  x = z^', 

x(l  - x)^  + ^ (1  - 9x)Dx  - | D = 0 (6.25) 

with  general  solution,  with  the  usual  notation  for  the  hyper- 
geometric  functions, 


D(x)  “ a-,F(l,  |,  x)  + a?x1/2F(|,  3’  I’ 

For  the  given  values,  the  two  independent  solutions  reduce  to 
the  finite  forms 


F(|,  3’  x)  = (1  - x) 


-3 


F(l,  x)  = - l(x2  - '.x  - 3)(1  - x)  3 


r 


(6.26) 


Therefore,  in  terms  of  6,  where  x = cos  Q, 


D(9)  = a2Cos  9 cosec^©  + a^Ccos^©  - 6 cos2©  - 3)cosec^©  (6.27) 


The  first  independent  solution  of  (6,23)  js,  by  comparison  with 
(6.22),  just  a constant  multiple  of  A(©),  This  is  otherwise 
obvious  since  if  D(t)  = constant  x A(t)  then  equations  (6.13) « 
and  (6,1*+)  become  identical  and  reduce  to 

Dtt  + (Y  + 3 )DDfc  + (y  + 1)D3  = 0 

2 -3 

Knowing  that  one  solution  of  (6,24-)  is  simply  z(l  - z ) 
it  is  possible  to  evade  the  general  method  (which  leads  to  the 
hypergeomctri  c e qua t f on  (6,25)  and  the  reduction  then  no  c 6 s s Q.T  y 
to  deduce  the  finite  forms  (6,26)  from  the  original  infinite 
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series)  by  seeking  a solution 


D(z)  = w(z)z(l  - z2)  3 


where  «(z)  is  to  be  determined  by  direct  insertion  into  (6.24-), 
This  gives  for  w(z) 

oj  2 > 

zz  _ _ 2(1  + z ) 


z(l  - z-) 


then 


w(z)  = — (zlf  - 6zJ  - 3)  + Kl 


if  K,  Kf  are  constants.  Therefore  for  the  second  independent 
solution  of  (6.24)  we  can  take 


D(z)  = (zh  - 6z3  - 3)(i  - z2)"3 


which  reproduces  the  second  term  of  (6,27)  with  z, 
When  also  expressed  in  terms  of  z, 


= co 


s e = Sx. . 


r— 3 


A - az(l  - z) 
at  = - * z(3z^  - 10z2  + 15) 


= - 5(1  - z2)2 


i I 


(6.29) 
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The  Flow  in  the  Physical  Plane.  Attention  will  be  con- 
centrated on  the  functions  given  by 

A(t)  = a z(l  - z2)'3;  D(t)  = Pz(l  - z2f3 


where 


at  = - | (3Z4"  - lOz^  t 15) 


then  clearly  for  v(y;t)  in  (6,11)  we  can  take,  with  p = - ayQj 


v = az(y  - yQ)(l  - z2)  3 


u(t)  is  given  by  (6.12)  rewritten,  in  virtue  of  (6,29), 

£ (log  V = - f jj|  = 2z(l  - z2)'1 

therefore  on  integration  in  terms  of  parameters  (jj,,u  ), 


- uQ  = ajxzCz^  - 3) 


(6.30) 


(6.31) 


The  formulae  (6. 30),  and  (6.31)  represent  the  complete  parametric 
solution  to  the  motion;  z being  given  in  terms  of  t by  (6,29), 

To  obtain  the  instantaneous  streamlines  (Fig.  3)  defined  by 


dy  _ v 
dx  u 


use  (6.30),  and  (6.31)  then 


ai  = (y  - yo)gCz) 


(where  g(z)  = ^ • for  u --  0)#  By  integration 

H(3  - z2)(l  - 22)3  0 


the  streamlines  are  obtained 


_ (x-xQ)g(z) 


(6.32) 
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since  z is  a function  oi  time  only. 

To  obtain  the  particle  paths,  the  simultaneous  equations 


at  * »<*>»  H * v(y:tJ 


in  which  z(t)  is  known,  must  be  integrated* 
Write 


2 s"1 


= v(y;z)  ||  = - 5z(y  - yQ)(l  - z2) 

after  using  (6.30),  and  (6.29).  On  integration  this  yields 

, 2 2/5 

1 - 7,  = K(y  - y0)  , 


(6.33) 


K being  a constant. 
Also 


2\2 


at  = u(z)  at  = f [uo  +“t‘z<z2  - 3)1(1  - z2> 

so  that  x may  be  expressed  as  a polynomial  in  z,  which,  in  con- 
junction with  (6.33)}  will  provide  the  particle  paths  in  terms  of 
z as  current  parameter* 

If  uQ  = 0 then 


dx  cfrf/r,2  , w,  „2\2 

= - 5p-(z  - 3) U - z ) 


and 


ay  _ 


y - y. 


^ dx  (z2  - 3) (1  - z2)3 


Substituting  for  z2  from  (6.33)  gives 

,.3  dy  , ..-1/5  f o . w N2/5|”~ 

-m-  K di-  - (y  - y0>  ]_2  + K(y  - y0} 
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therefore 

- — ^-r  = \ (y  - y 0)6/^  + | (y  - y0)8/^  + K‘  (6.34) 

5p,K3  3 0 8 0 

the  constants  of  Integration  (K,  K' ) being  used  to  specify  the 
ccF  particles,  by,  for  example,  their  positions  at  some  specified 

K 

time.  We  will  consider  a < 0,  < 0,  On  putting  r)  = y - y0 

the  particle  paths  may  be  written  parametrically  in  terms  of  t) 
to  rationalize  (6 *34) , 

- ilg.  . - n6(8  + 3Kr)2)  + K* 

5}iK3 

where  u is  a constant  for  the  flow  as  a whole  but  where  K,  K1 
are  used  to  parameterize  the  od  particles.  The  fluid  particle 
(K,K'  ) is  situated  at  [x(t}),  y(  rj)3  at  the  time  t(r))  given  by 

I 2 2 If  2 

3 at  = - yi  - KT)  • (3K  n +^Kn  8) 

(obtained  by  eliminating  z between  (6,29)  and  (6.33).  In  the 
present  solution  |z  | 1 therefore,  since  1 - z = Kt)  ,the  para- 

meter K must  be  positive  for  all  particles.  It  follows  that  the 
current  parameter  rj  is  confined  to  the  range  (-  K^2,  +K^'//2)» 

By  first  drawing  the  separate  curves  of  (x,y)  against  rj  and  after- 
wards combining  them  to  obtain  the  particle  paths  given  by  the 

-1/2 

x against  y curves.  At  zero  time  r)  = K and  z = 0 therefore 
u = v = 0 everywhere.  Thus  the  fluid  is  instantaneously  at  rest 
at  zero  time  when  the  particles  begin  to  advance  towards  the 
line  y = yQ,  arriving  there  at  t = -8/3 a ( > 0)  when  n = 0,  It 
is  easily  verified  that 
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3u  _ 3n°  av  , „ 6v  _ a 2 

et  " — ’ H + * 3y  “ ~ 


5Kn  " 5KWV 

therefore  the  acceleration  becomes  infinite  when  r\  - 0.  The 
curves  of  x and  y against  t)  and  the  final  (x,y)  curve  for  this 
solution  are  denoted  by  the  full  lines  in  Figures  12,  13  and  14* 
The  present  analysis  cannot  yield  a solution  for 
t > - 8/3  a or  for  y < yQ. 

In  the  equation  preceding  (6.20)  it  would  be  equally 
legitimate  to  use 


= ,u2or2/5  + 

and  the  relevant  equation-  now  become,  if  a = p\|> 
= $ 0(z2  - 1)^  2 
A - cz(z  -1) 


-6/5 


3at  = zOz1*  - 10z2  + 15) 


af|  = ^ 

v(yst)  = (y  - y„)A(t)  =az(y  - y JCz^  - 1)  3 

V o 

u(t)  = u_  + ap,  s(z2  - 3) 

The  instantaneous  streamlines  are  now  given  by,  when  \iQ  = 0, 

+(x-x  )g(z) 

y - y0  = b 

|!g(z)  = (z2  - 3)  ■L(z2  - 1)"^ 


and  the  particle  paths  by 
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y - yD  = n z2  - i = kv 

= ri6(3Kr,2-  8)  + K» 

fjj,  KT 

3 at  .[3^-^+  8] 

where  [ x(  t})  , y(r])3,  z ( rj)  give  the  position  ana  speed  as  above  at 
time  t ( rj)  of  the  particle  (K,  K*)»  We  may  choose  yQ  = 0 without 
loss  of  generality. 

This  second  solution  applies  when  j z | 2.  1 anc*  therefore 

the  parameter  K must  be  positive  for  all  particles  since 
2 2 

z - 1 = Xr,  . Unlike  the  first  solution,  however,  the  parameter 
may  take  on  all  real  values.  In  the  diagrams  (Figs.  12,  13, 
and  14)  the  first  solution  is  shown  by  the  full  lines  and  the 
second  solution  by  the  broken  lines.  The  first  solution  ceases 

Q 

to  apply  when  the  time  t = | — | is  reached,  and  the  second 

~>a  a 

solution  only  applies  for  t > j it-j  . Therefore  a solution  valid 

3c 

for  C <,  c <oo  may  be  obtained  by  concentrating  attention  on  the 
section  0A(  Y]  >_  0)  of  the  first  solution  in  conjunction  with  sec- 
tion OB  (n  <0)  of  the  second.  Initially  only  the  upper  half 

plane  is  occupied  by  fluid  instantaneously  at  rest,  corresponding 

—1  /2 

to  point  A (z  = t = 0,  n = K ).  When  point  0 (n  = 0)  is 

reached  the  fluid  is  moving  with  velocity  approaching  oo  and  the 

acceleration  is  infinite  at  time  | JL|»  If  the  second  solution 

pa 

OB  is  then  considered,  the  fluid  concentrated  along  the  x-axis 
expands  into  the  lower  half  plane  and  eventually  recedes  to 
infinity  parallel  to  this  axis.  The  diagrams  show  only  the  motion 
of  a typical  value  of  K and  when  K>  = 0,  The  constant  K1  merely 
serves  to  provide  the  x-wise  position  of  the  fluid  elements  at, 
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say,  t = 0. 

The  phenomenon  which  takes  place  on  the  axis  is  not  a 
simple  shock.  In  the  initial  motion  the  advancing  layers  of  f3.uid 
overtake  themselves  and  cause  a concentration  of  fluid  moving  at 
high  speed  across  the  axis,  thereafter  reducing  speed  during  the 
subsequent  motion  in  the  lower  half  plane,  but  accelerating, 
again,  to  the  r.Lght  as  t —*>00, 

Conclusion.  The  more  important  results  of  this  work  may 
be  summarized  as  follows • 

1.  In  unsteady  potential  flow  of  the  simple  wave  pattern, 
the  lines  of  constant  (q,0)  arc  straight.  However,  in 
the  general  flow  q(X;t),  the  pressure  and  density  and 
acoustic  speed  are  not  constant  along  such  lines,  except 
in  the  simplified  case  a(X).  In  all  cases  these  lines 
may  constitute  a system  of  parallel  lines  parallel  to 

a fixed  direction  but  may  never  be  concurrent.  In 
gonoral  the  X-lines  envelope  some  curve,  which  grows 
linearly  with  time  in  the  simplified  wave, 

2,  In  the  general  wave . the  solution  for  © and  $ may  be 
expressed  in  terms  of  four  functions  of  the  parameter 
X(q;t)j  these  functions  being  chosen  so  that  together 
they  satisfy  a Monge-Arnpere  equation  for  X(q;t).  This 
representation  is  sufficient  to  show  that  the  X-lines 
can  never  be  concurrent  in  the  physical  plane  and  that 
it  is  net  possible  to  obtain  any  flows  of  an  oscillatory 
nature. 
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These  deductions  are  in  contrast  to  the  steady  flow 
phenomena.  For  steady  potential  flow  the  \-lines  may  be  parallel } 
concurrent  or  possess  a non-degenerate  envelope.  If  the  steady 
iso-energetic  flows  become  rotational  then  the  lines  never  possess 
a non-degenerate  envelope  and  must  either  be  parallel  or  con- 
current, In  contrast  the  latter  alternative  is  the  very  case 
which  is  excluded  when  the  simple  wave  becomes  unsteady.  It  is 
useless  therefore  to  attempt  to  obtain  an  unsteady  simple  centred 
wave,  just  as  it  is  impossible  to  generalize  the  steady  rotational 
simple  centred  wave- 
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